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ABSTRACT 


Initial yield surfaces of unidirectional fiber reinforced 
composites subjected to combined longitudinal, transverse 
normal, and longitudinal shear loads are calculated. The 
composite is composed of filaments arranged in a periodic 
Square array embedded in a matrix material. The constituent 
Materials are assumed to be homogeneous and isotropic, while 
the composite is assumed to be macroscopically homogeneous 
and transversely isotropic. The finite element method, 
using linear strain triangles, is employed to calculate the 
stresses throughout the composite. Von Mises yield criterion 
is used to calculate the elastic limit of the local micro- 
scopic combined stresses in the composite. A parametric 
evaluation is carried out by changing individual constituent 
properties and evaluating the effect of this variation on 
the composite properties. Yield surfaces for several 
functional composites are included, as well as a computer 
program, in FORTRAN IV language, for the calculation of a 


yield surface. 
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I. INTRODUCTION 


A. COMPOSITE MATERIALS 

The requirements of modern technology, especially 
that of the aerospace industry, have necessitated the 
development of materials possessing much higher strength 
to weight ratios than common structural materials. This 
result has been achieved through the introduction of com- 
posite materials, that is two or more constituent materials 
combined in such a way that the gross properties of the 
composite are favorable combinations of those of the 
constituents. This indicates that the desirable properties 
such as strength/weight and stiffness/weight ratios are 
increased while less desirable properties such as .low 
ductility are minimized. 

A specific type of a two constituent composite is a 
unidirectional fiber reinforced composite. This consists 
of continuous lengths of a stiff reinforcing fiber imbedded 
in a ductile matrix material. All the reinforcing fibers 
are Oriented in the same direction within Single layer of 
material. Generally several layers, or laminae, are bonded 
one on top of another to form a high strength laminated 
composite. 

Each of the lamina may be considered as an anisotropic, 
homogeneous material. The transverse properties, i.e., 
those in a direction perpendicular to the filament orienta- 


tion, are weaker than those in the longitudinal (fiber) 
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direction for a single lamina. To counteract this weakness 
in a laminated composite, adjacent lamina are arranged such 
that they have different filament orientation. The funda- 
mental unit is therefore the single layer of unidirectionally 
reinforced materials. It is this type of basic layer that 

is analyzed in this investigation. 

The physical dimensions of such a lamina can vary 
considerably, although the largest is generally but a fraction 
of an inch thick. The fibers vary in diameter from .0001" 
mor glass and .0003" for graphite to .004" for boron. Lamina 
thickness depends on the materials used and the intended 
purpose. An individual lamina can be produced in nearly any 
desired thickness. However standard laminae are becoming 
popular. For example, most graphite comes in a standard .0l" 
thickness. Allowing for the presence of the matrix material 
such a lamina may contain in excess of 25 fibers across its 


thickness. 


B. DEFINITIONS 

The terminology of composite materials is often confusing 
due to a lack of precision in defining specific terms. For 
example,at times it is unclear whether a term such as tensile 
strength refers to the ultimate strength of a material, i.e., 
its failure point, or to the point of initial yield. In 
composite materials it is necessary to be especially clear 
by noting whether a term refers to a composite property or 
to a property of a constituent material. In order to avoid 


such difficulties the definitions of this section will be 


adhered to in this study. 
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1. Macroscopic Quantities 
These are quantities of the composite material. 
They are calculated by assuming that the composite consists 
of anisotropic, homogeneous layers. 
2. Microscopic Quantities 
These quantities are determined on a point by 
point basis throughout the composite. The different 
mechanical properties of the matrix and filament materials 
are considered in the calculation of these quantities. To 
see the relationship between macroscopic and microscopic 
quantities consider Fig. (1). Figure (1a) represents a 
body with microstress, Ore In Fig. (1b) the macrostress 
S in the x-direction obtained from the microstress field 
of Fig. (1a) is 


5 = 2 || c. (y) dy (el) 
x a 


O 
E. is the macrostress for the face x = b, and represents 


the average value of the microstresses. 
3. Initial Yield Strength 
That combination of macroscopic stresses at which 
a point in the composite first reaches its elastic limit, 
i.e., no plastic strain has occurred in the composite. 
This study is restricted to plane macrostress 


fields with components Ва апа 5х2, Where the z axis 


z! 
is along the length of the fiber. Specific cases of 
initial yield strength associated with this investigation 


are 
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a) Microstress Field c (b,y) 
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5 
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x 
О b 


b) Macrostress 5. 
<` 


Figure (1). Relationship of Microstresses and Macrostresses 


a) Longitudinal Yield Strength, when the only 
macroscopic stress is in the axial direction, i.e., 
S, A 0, S. = SE 0. 


b) Transverse Yield Strength, when the only macro- 


Scopic stress is in the transverse normal direction, i.e., 
S. x 0, S = S = 0. 

c) Longitudinal Shear Yield Strength, when the 
only macroscopic stress is in the longitudinal shear 


earreetion, i.e., SB * 0, 5. = 5, = 0. 
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4. Second Deviated Stress Invariant, Ј, 
The second deviated stress invariant at a point 


is defined as [1] 


"nm | (1-2) 


3 2 
(cy ¬ Og) + (o 
This quantity can be related to the yield strength of a 
material. 
5. Effective Microstress, O 


The effective microstress, 0, is proportional to 


the square root of the second deviated stress invariant, [1] 


zu d _ 2 _ 2 _ 2 
O r4 [< су) + (o; с.) + (97 ο) 
+ 6 cy + ΘΙ + 2 3 (1-3) 


6. Specific Strength 
The ratio of a particular initial strength term, 
as defined above, to the specific density of the composite. 
The strength considered must be specified, for example, 
specific longitudinal strength. This may be applied to 
composite materials as well as single phase materials. 
7. Specific Modulus 
The ratio of the Young's modulus of a material 
to the specific density of the material. In the case of 
anisotropic materials it must be specified which modulus 


is being considered, such as the specific transverse modulus. 
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8. Ultimate Strength 
That point in macroscopic stress space at which 
either a matrix element or a filament element actually 
fails, i.e., breaks,cracks, or separates. Due to the 
brittle nature of filament materials the initial yield 
strength and ultimate strength usually coincide for filament 
Materials. 
9. Initial Yield Surface 
In a general three-dimensional system the stress 
state is defined by six stress components, three normal 


stresses, S т апа S,,and three shearing stresses, S 


xy' 
Thus an initial yield surface may be repre- 


xt 


S and S 


pez" y z ° 


sented in macrostress space. Any macroscopic stress state 
within the surface means that the effective microstress at 
every point in the composite is less than that required for 
yleld. Any stress state lying on the surface means that 
some point in the material has an effective microstress 
equal to the elastic limit. For brittle materials without 
a plastic region loading states outside this surface are not 
possible. For ductile materials stress states outside the 
initial yield surface (within the limits of the ultimate 
strength of the material) will result in plastic deformation. 
The three-dimensional yield surface considered in 
this study is a specific case of the general initial yleld 
and 


surface described above. Macroscopic stresses, S S 


у! xy’ 
Syz are set equal to zero in the calculation of the yield 


surface. The resultant surface represents the yield surface 


due to plane macrostress loading conditions. 
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10. Stress Concentration Factor 
A stress concentration factor, denoted as SCF, 
is generally defined as the maximum microstress in the 
composite to the applied macrostress. In this study three 
particular cases of stress concentration factors are 
investigated. 


a) SCFX. Transverse loading stress concentration 


factor, associated with macrostress state (S,, S,, 5,,) = 
(S., 0, 0). 
(9х) 
ЅСЕХ = E (I-4) 
X 


b) SCFZ. Longitudinal loading stress concentration 


factor, associated with macrostress state (S,, S,, Syz) = 
pos, 0). 
(55) 
max 
SCF = m TE) 
2 


c) SCFXZ. Longitudinal shear loading stress 


concentration factor, associated with macrostress state 


са = Ru. (1-6) 

11. Transverse Stiffness, Em 
The transverse stiffness of a composite is the 

macroscopic modulus of elasticity in a direction perpen- 


dicular to the fiber direction. 
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12. Longitudinal Stiffness, 
The longitudinal stiffness of a composite is the 
macroscopic modulus of elasticity in the fiber direction. 


13. Longitudinal Poisson's Ratio, ντ 


The longitudinal Poisson's ratio is the measure of 


the contraction in the transverse direction due to tension 


. . : : € 
ο. ὃς longitudinal direction, 1.е., ντ = = , for the case 
Ζ 


14. Transverse Poisson's Ratio, V 
—— = TA 
The transverse Poisson's ratio is the measure of 
the contraction in the transverse direction due to tension 


in the same plane, that is v, = Ey Ey for the case 


T 
(S 


о. 0, О. 


Y Z 
15. Shear Modulus, G 


x 


The composite shear modulus is the modulus of 
elasticity in shear due to longitudinal shear loading. As 
the transverse composite shear modulus is a dependent 
quantity any reference to composite shear modulus in this 
study will mean the longitudinal shear modulus unless 
specifically stated to the contrary. 

16. Filament Volume Fraction, Ve 


The filament volume fraction is the ratio of 


filament material volume to total material volume, 


το Filament Volume 
f Total Volume (I-7) 


Note that 0 < νε ο. 
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l7. Matrix Volume Fraction, v 
———v— nŠIO 
The matrix volume fraction is the ratio of the 


matrix material volume to the total material volume, 


у = Matrix Volume 
m Total Volume (1-8) 


Again, 0 < v, « 1.0. For a material without voids the sum 
of the matrix and filament volume fractions must equal 
unity, 


C. BACKGROUND 

When a tensile test is performed on a composite material 
specimen to determine an ultimate tensile strength the 
value obtained is a macroscopic property. When design 
specifications require certain properties these also are 
macroscopic properties. To accurately determine the behavior 
of a composite and to bridge the gap between design require- 
ments and design practice, macroscopic behavior must be 
related to microscopic behavior within the composite. 

Numerous analytical approaches for the determination of 
composite properties have been utilized. These include 
strength-of-materials, stress function, finite difference, 
and finite element techniques, among others.  Chamis and 
Sendecky) [2] present a comprehensive review of the various 
analytical methods. 

Specific examples of micro-stress analysis include the 


longitudinal shear loading solutions of Adams and Doner, [3] 
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and Tsai, Adams and Doner [4]. The longitudinal load 
case was considered by Bloom and Wilson [5] using a series 
solution technique. Transverse normal loading was investi- 
gated by Tsai, Adams and Doner [4] , and Adams and Doner [6] 
using finite difference analyses, and by Foye [7] using a 
finite element method. Lin, Salinas and Ito [8] етр1оуеа 
finite element techniques to consider the problem of 
combined loading, by superposing the effects of longitudinal, 
transverse, and longitudinal shear loading. Lin et al 
calculated an initial yield surface based on this combined 
loading analysis and later extended their work into the 
pBESStic region, [9]. [10]. Adams [1] also considered inelastic 
loading but restricted his study to transverse normal loading. 
The combined loading problem has received the least 
attention though it is the most practical from a design stand- 
point. It is of great importance to the designer to know 
what combination of macrostresses will result in a micro- 
stress distribution sufficient to cause yield. The initial 


yield surface provides this information. 


D. SCOPE AND OBJECTIVES 

The analytical approach used in this investigation was 
originally developed and programmed by Lin, Salinas, and 
Ito. [8] The author's objectives have been 

1) Modify the program to allow use of an overlay 
system structure in order to reduce the required computer 
storage. 

ii) Substitute the use of random access disk unit for 
the original tape storage methods. 
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iii) Perform a parametric study of the effects of 
constituent properties on the properties of the composite. 
iv) Compare several realistic material combinations to 
evaluate their structural efficiency. 
v) Provide selected initial yield surface plots for 
design assistance. 
vi) Make some general observations relating composite 


behavior to constituent materials. 
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II. METHOD OF ANALYSIS 


The analytical method described in this section was 
developed by Lin, Salinas, and Ito [8]. The present discus- 
sion is a general presentation of the basic principles of 
that analysis. Appendix A is a detailed description of 
the method. 

The basic assumptions necessary to the analysis are, 

i) The filament and matrix materials are isotropic 


and homogeneous. 


ii) The filament and matrix materials are perfectly 
bonded 9 4 
It should also be noted that the macroscopic stress 


state is one of plane stress, with components Sy S and 


2? 
ES 

The microscopic stress distribution within a unidirec- 
tionally reinforced composite depends upon the arrangement 
of the filament material within the matrix. In order to 
facilitate this EE the filament array is assumed to 
be a doubly periodic rectangular array. In such an array 
the two-dimensional (x-y) surface may be divided into 
identical segments called basic blocks, of width 2a and 
height 2b as shown in Fig. (2). Away from the edges of the 
composite the microscopic stress distribution should be 


identical within each block. Due to symmetry within a basic 


block only one quadrant need be analyzed, as shown in Fig. (3). 
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Eure (2). Rectangular Array of Filaments and a Basic Block [5] 


Y 


Houses). A Quadrant or the Basic Block [8] 
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In order that no voids be formed and that adjacent 
blocks neither separate or overlap it is necessary that 
the rectangular basic blocks remain rectangular during 
deformation. In order to produce a combined macroscopic 
plane stress state, as well as maintaining the rectangular 
form of the blocks, four specific load cases are defined. 
These four cases, shown in Fig. (4), are superposed to 
give plane stress and boundary compatibility. 

Load cases I and II are plane strain states with 
ES = 0. Load case III is a generalized plane strain case 
with unit strain. Load cases I, II, and III are superposed 
to give S, = 0. Load case IV is a state of longitudinal 
shear. All four load cases have Sxy = Syz = 0. The boundary 
conditions for each problem are shown in Fig. (4). 

These four load cases are solved by the finite element 
method using linear strain triangles. Felippa 19] presents 
a detailed discussion of finite element theory and a 
derivation of the linear strain triangle. The finite 
element model utilized for the circular filament cross 
section problems is shown in Fig. (5). This model was 
adapted from that of Lin et al [9]. In this work convergence 
tests were conducted to establish the validity of the model. 

Each load case yields a microstress distribution and 
associated macrostresses. The microstresses are designated 
by В. with subscripts determining the stress component 


І I 
and superscripts denoting the load case. For example ij 


is the microstress 0;. at a specific point associated with 


j 
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problem 1. Macroscopic stresses are represented by S^ 
s: 
with subscripts and superscripts similar to T 
Appendix A treats the solution of the individual 
problems and demonstrates how the results may be combined 


to relate the combined loading microstresses at a ρου EO 


macrostresses as 


[о] = [a] [5] (11-1) 
See Appendix A, Eq. (A-13), for the explicit form of these 
matrices. 

As noted, the microstresses vary from point to point. 


Therefore the A matrix, evaluated as 


[a] = [6] [в] (11-2) 
is also point dependent. 

In this thesis the Von Mises yield criteria was taken 
to govern initial yield of the composite. In terms оГ 
microscopic stresses, yielding occurs when [1] 

J = cae - our t o L |o +o + و رت‎ - Oy 9z - 0) 


1 2 
= 4 II-3) 
5 ( 


where Yo is the elastic limit of the material for the tensile 
stress state. 

Using Eq. (II-2) to relate microscopic and macroscopic 
stresses Appendix A shows that the yield equation, Eq. (Е, 


may be written as 


28 





2 
E S, + C 


1 4" xz O 


2 о E 2 2 
29x89, + Сз5х + C,S = S Y (II-4) 


The terms C Cor C4, and σι are second order polynomials 


1? 
of the components of the A matrix and are functions of 
position. For given values of Sy. and S,zı Eq. (11-4) пау 

be solved for the value of S, necessary to cause yield. 

Each point within the fundamental block will yield a dif- 
ferent value of S Each combination of Sy and S,, yields 
two values of S, at each point, a maximum and a minimum. 

The minimum of the maximum values and the maximum of the 
minimum values determine two yield points for the given 
values of ER and So" Continuing in this manner the initial 


yield surface for a composite subjected to a macroscopic 


plane stress state is calculated. 
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ТТТ. EVALUATION OF NUMERICAL RESULTS 


А. INITIAL YIELD SURFACE CHARACTERISTICS 

The two-dimensional plots of the yield loci of this 
investigation have the general form of Fig. (6). Each 
locus represents a specific value of longitudinal shear 


stress, S with the outermost locus corresponding to 


τ. 


E - 0, as seen in Fig. (6). 

Each point on a yield locus indicates that some point 
within the material would yield under the given macrostress 
state. The analysis used in this investigation provides 
the location of the area of yield associated with each 
point of a yield locus. These area locations refer to the 
triangular elements of the finite element grid, Fig. (5). 
Figure (6) is the yield surface for a 50$ Boron/6061 
Aluminum composite. The elemental areas that first reach 
yield under specific loading conditions are listed in Table 
I. The corresponding macrostress states are noted on Fig. 
(6). Proceeding in this manner a particular triangular 
element could be associated with every point on each yield 
locus. This information could be useful to a designer in 
determining if failure is associated with a particular 
region. 

On initial consideration it seems somewhat unusual 


that a yield locus is not symmetric about the S, and S axes. 
2 


*. 


To explain this behavior first recall the yield equation of 


the previous section, Eq. (11-4). 
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The exact definition of the terms Ci: Car Car and Ca can be 


ascertained from Eq. (A-22). 





TABLE I. INITIAL YIELD FOR 5» = 0 
Point S (ks1) S (ksi) Yield Location 
E 2 (Element) 
A 0 114 66 
А' 0 -114 66 
Β вето ΙΗ; 67 
B’ TOSS US 88 
С 2030 +100 67 
ον 30 +.45 88 


As stated, the points of the initial yield surface are 


calculated by setting S,, and S, to specified values and 


calculating the value of 5, resulting in a microscopic stress 
sufficient to cause yield. Assume that the macrostress 
S 


combination (S x S21) is known to be such a yield point. 


xc 
It will be shown that for different combinations of E and 
521, Terr G S 5,1) etc., a specific type of symmetry 


exists. 


Bor (Sy; 5,1) Eq. (III-1) becomes 


y 2 (III-2) 


2 2 2. 
С1521 + С25215х + C38y ^ + CgSgzó= = Y, 


w |} 
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E = 5 |) Eg. (III-1) is 


2 2 2m 1 2 
Ci9z1 T C2535315, + C35, + Ca?xz = = Yo (III-3) 


Equations (III-2) and (III-3) are identical and show 
that yield occurs at the same point and for the same 
absolute values of macrostresses in quadrants I and III 
of the yield locus, as shown in Fig. (7). This indicates 
symmetry about the Бс axis, perpendicular to the (S. ,S;) 
plane. 


For combinations of (Sy r 5.1) such that their product is 


negative, (S,, E or (-5.., 5,1), Eq. (III-1) becomes 


2 _ 2 
С 521 ео 


1 
S 2= _y2? (III-4) 
Ax 3 Q 


Ewbaring Egs. (III-2) and (III-4) it is seen that if (Szr 


S I IS a valid solution of Eg. (III-I) as was assumed, 


x7 


then (S,z, 5 тг) and (S -S 521) cannot be a solution. 


κ’ xz” x 
For a positive S, and писал so Solution ОЕ БАТЕ) 
will yield a different minimum S,, say $5,2. This corresponds 
to a different point in the basic block than does the 
solution for quadrants I and III. The combinations (S,, -5,„5, 
Syz) and (-S,, S22, =) identically satisfy Eq. (III-1) 
establishing the equality of yield points in quadrants II 
and IV of Fig. (7). The symmetry of the yield locus about 
the S,7 axis is thus established. 

Another interesting phenomenon apparent in some yield 


surfaces is the occurrence of square ends on the yield loci 


eorresponding to small values of S,,. This is shown in Fig. 
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(8). Generally this indicates that yield is occurring in 
the filament material rather than the matrix. Note that 
the S, scale e Eon οσοι larger than the S. scale. 


The longitudinal stress values, S in the vicinity of the 


7!’ 
square ends are many times larger than the transverse, S,, 
Memes, The S, values are so dominant that increasing the 
transverse stress from point A to point B of Fig. (8) has 


little effect on the yield surface. Once the macrostress 


state B is reached yield occurs in the matrix material. 





Figure (7). Symmetric Nature of a Yield Locus 


As stated the yield locus plot is a two-dimensional 
representation of a three-dimensional surface. Thus the 
Separation between successive yield loci is an indication 
of how steep the slope of the yield surface is in that area. 
The dark areas of the plots indicate an extremely steep slope, 
such that adjacent loci actually plot very close to one 


another. 
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B. EVALUATION OF COMPOSITE MECHANICAL PROPERTIES 
1. Composite Strength 

Composite strength properties were determined 
directly from the yield surface associated with each 
particular composite. 

i) Longitudinal yield strength was taken as 
points S, on the initial yield surface when S, „= Sy= 0. 

ii) Transverse yield strength was taken as points 
on the initial yield surface when S; = Sy, = 0. 

iii) Longitudinal shear yield strength was taken as 
the point S,,, on the initial yield surface when Sy= S, = 0. 
2. Composite Elastic Constants 

Composite elastic constants were calculated from the 
results of the individual longitudinal, transverse, and 
longitudinal shear loading problems. The macroscopic 
stresses associated with each loading case were used in 
evaluating the transversely isotropic stress-strain relations 
developed in Appendix B. As discussed in this appendix a 
transversely isotropic material is an orthotropic material 
with a single plane of isotropy. It should be noted that 
the assumption of transverse isotropy limits the applicability 
of this method of determining the elastic constants of a 
composite. Composites are not generally transversely 
isotropic. The particular case treated in this investigation 
is that of a square array of filaments. For this particular 
fiber distribution the composite properties approximate those 


of a transversely isotropic material. For more general 
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filament arrangements, such as a rectangular array, further 
consideration must be given to the applicability of these 
stress-strain relations. The analytical method employed 
to calculate the displacements, stresses — strains is 
not effected by this restriction however. The calculation 
of these quantities is applicable to a general rectangular 
filament array. 

For convenience the transversely isotropic stress-strain 


relations are repeated here. 


l y 
a) £ — = mE - 1, 
x (S VS) — с 
En x T y Er, Ζ 
ν 
ο ος τς νο) - b 
Y S 
ν 
E IA 1 E 
ο) ES É Tm * Sy) + y S; (III-5) 
L i 
1 
а € ME 
E =1 ۽‎ 
y2 7 
σι vA 
1 


In the transverse loading case, problem 1, the macro- 


scopic strain, є corresponding to the unit displacement 


x 


imposed, was calculated. The macroscopic strains Ey and 


є, аге equal to zero. 
In the longitudinal loading case, problem 111, the macro- 


scopic strain εσ is equal to unity. In this problem 
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Equations (III-5), evaluated for these two loading 
conditions, were solved for longitudinal and transverse 
E ν 


composite stiffness and Poisson's ratios, E ν 


ΤΊ Τ᾽ 
The longitudinal Poisson's ratio was evaluated directly 


from Eq. (III-5c) in terms of load case I as 
πο === (III-6) 


Longitudinal stiffness was evaluated from Eq. (III-5c) in 


terms of load case III as 


III III 
E + gill 


τ 5 LU. y ) (III-7) 


Equation (III-5b) may be evaluated in terms of load cases 


I and III and the resulting eguations solved simultaneously 


give Д 
III 
e Ey, Sy sl. TIL sI 
= E Y ὃν x £ 
I M iS ЕТ eo) 
L Smash. 5 
7. x 7. x 
I 
E, EE sini 
A MEE n М» |ы .ϱ (III-9) 
T E TI = 
т SET ST I I 


The longitudinal composite shear modulus, G, was obtained 
from the results of the longitudinal shear loading case, 
problem IV. The shear strain corresponding to the imposed 
unit displacement was calculated and used in conjunction 


with the longitudinal shear strain Eq. (III-5f). 
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3. Specific Strength and Specific Stiffness 
Specific strength and specific stiffness were 
calculated by dividing the respective properties of the 
composite by the composite density. Densities of the 
various material combinations were calculated from the 
specific weight data of Table II, which lists properties 
for all materials used in this study. 
4. Stress Concentration Factor 
Three specific values of the stress concentration 
factor are calculated in this investigation, corresponding 
to three specific loading cases. All comply with the 
general definition of the ratio of the maximum microscopic 
stress to the macroscopic stress. 


i) SCFX refers to the transverse loading problem, 


problem I. It is defined as the maximum microscopic stress, 
бу, divided by the average surface traction, 5х. From Eq. 
EU (ох) 

SCFX = Tc (III-10) 


IP SE rerers CG Ehe longitudinal loading problem; 
problem III. It is defined as the maximum microscopic 
stress, O,, divided by the average surface traction. From 
Ед. (1-5), 


(σ 
= max 
5СЕ2 = TS (IIT) 


111) SCFXZ refers to the longitudinal shear loading 
problem, problem IV. It is defined as the maximum micro- 
scopic shear stress Ох, divided by St the average shear 


stress of this problem. 
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С. COMPARISON OF NUMERICAL RESULTS 


In order to establish the validity of the analysis 
used in this investigation comparisons were made with avail- 
able theoretical and experimental results. It is not the 
intention here to establish the superiority of one method 
of analysis over another but rather to ensure that the 
solutions are compatible. 

1. Transverse Stiffness, En 

Using a finite difference approach Tsai et al [4 | 
determined the transverse stiffness of a 40% fiber volume 
Boron/6061 Aluminum composite as En = 20 x lo psi. mhe 
present analysis gives E, = 20.97 x 10° pasi, 

Adams and Doner, [6 | also using a finite difference 
analysis, calculated a theoretical value of transverse 
Stiffness for a 47% E-glass/Epoxy model as 1.55 x 10° psi, 
and reported an experimental value of 1.30 x 10° psi. This 
analysis gives 1.42 x 10° psi. Adams and Doner also list 
a comparison between their results for Glass/Epoxy and 
several other analyses, including finite element, complex 
variables, and other solution techniques. These results 
are compared to the present analyses in Table III. 

Transverse stiffness for Boron/Epoxy and Glass/Epoxy 
also fell within the variation bounds calculated by Dow 


et al πα. Further comparisons were made with results of 
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TABLE ТТ. MATERIAL PROPERTIES 


Specific Wt. 


_ Material! E(psi) ν Y (psi)? αρ) __ 
1. Boron 60x10° 0.2 4.8x10° . 085 
2. S-glass 12x10° 082 5.0x10° .090 
3. Е-91аѕѕ πο ο NOS ουν» 4.0x10° ‚ 090 
4.  MODMOR I 60x10° 0.25 OR 2072 

(graphite) 
5.  MODMOR II 40x10° 0.25 ZONE .063 
(graphite) 
6. THORNEL 25 25x10® 0.25 1.8x10° .052 
(graphite) 
7.  THORNEL 40 40x10° 0.25 2.5x10° . 052 
(graphite) 
8. 6061 Aluminum 10x10® 0.3 .346x10? .100 
9. 2024 Aluminum πο 034 130105 .100 
10. NARMCO 2387 .5x105 0.3 Π.Ο: 72045 
(epoxy) 
11. 4617 Epoxy .78x10* 0.34 .800x10° .045 


1. Materials 1 througn 7 are filament materials. 
Materials 8 through 11] are matrix materials. 
25 Y. is the elastic limit of the material for the 


tensile stress state. 
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Ekvall 16], Whitney 17] and others. The majority of 
these results are in the form of plotted curves which 
makes direct numerical comparisons difficult but does 


provide a general comparison of results. 


TABLE III. COMPARISON OF TRANSVERSE NORMAL LOADING SOLUTIONS 


E-/E=20 νε-.2 %д=.34 ν.-.55 
Analysis Transverse Stiffness, Em 
Adams and Doner [6] 1.86x10° psi 
Wilson [13] 1.86x10% psi 
Pickett [14] 1.85x10% psi 
Foye [7] 5510551 
Present analysis (473) 1.42x10% psi 
Present analysis (60%) 2 MOS 
Present analysis (50% estimate) | 61.10 psi 


2. Composite Shear Modulus 
Adams and Doner Bis using a finite difference method, 
have calculated numerous theoretical values of composite 
shear mođuli for specific material combinations. For 
Glass/Epoxy they obtained a shear modulus of .417x10° psi 
and .589x10% psi for 40% and 55% models respectively. The 
present analysis yields .447x10° psi for a 47% model of the 


same material. 
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For Boron/Epoxy Adams and Doner obtained .441x10*? psi 
and .658x10% psi for 40$ and 50$ models, compared to the 
present values of .476x10% psi for a 47% model. 

Direct numerical comparisons were available from 
Blakslee et al fi 8] who utilized Whitney's micromechanics 
model 17 based on an Airy stress function. The particular 
graphite investigated was "Thornel" (a commerical graphite 
produced by Union Carbide). For Thornel 25 in a 50% model 
Blakslee reported a theoretical value of the composite shear 
modulus of .518x10° psi and an experimental value of .60x10° 
psi. This compares to a present value of .531x10° psi. 
Similar values for a 60% model of Thornel 40 are .931x10° psi 
theoretical and .740x10° psi experimental, compared to a 
Present value of .768x10° psi. These comparative values for 
Thornel are summarized in Table IV. 

Graphical comparisons were available with Ekvall [16], 


Tsai et al, [4] and Whitney |7] . 


TABLE IV. SHEAR MODULUS COMPARISON FOR THORNEL/EPOXY 


Shear Modulus, G(psi x 10 9) 


Analysis THORNERT 25? THORNEL 40** 
Blakslee [18] .518 .931 
Experimental [18] „600 „740 
Present Analysis JD .768 


* THORNEL 25 E,/Eg- 32 v7.25  wvQ7.34  vg-.5 


** THORNEL 40 Eg/E - 51  vg-.25 v7.34 νε-.65 
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The graphical results of Fig. (9) were compiled 
by Reidinger et al [19] . This compares the values derived 
by Tsai [20], Dow et al [21] , and Ekvall [16] for Glass/Epoxy. 
The results of this analysis are marked by X. Both Tsai and 
Dow et al obtained their results from consideration of 
variational bounds. Tsai however employs a "contiguity" 
factor dependent on the amount of contact between filaments. 
The present analysis and that of Dow et al assume that the 
filament is completely surrounded by matrix material. This 


accounts for the closer correspondence of the present results 


to that of Dow et al. 





220 
—— Fsal [20] | 
— Dow et al [15] : 
z l; 
A 1.5 ЕН 0.5] 2 а 
3 x Present Analysis [и 
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E ] 
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> τς 
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ο ο NES 1,5 .06 7. „8 


Fiber Volume Fraction (ve) 


E.-10.6x10* psi E,-0.5x10* psi vg-.22 vg7.35 


Figure (9). Shear Modulus, G, vs. Fiber Volume Fraction,v. 
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3. Longitudinal Stiffness, Er 
Longitudinal loading has not been considered as 
frequently as transverse or shear loading, resulting in 
fewer numerical comparisons. Blakslee [18] again gives 
theoretical and experimental values for Thornel/Epoxy 
composites. For the 50% volume Thornel 25, he reports a 
Eheoretical value of 12.78 x 10% psi, and an experimental 
meme of 11.07 x 10% psi. This compares to a value of 
11.96 x 10% psi in this study. For 65% Thornel 40 similar 
Besults are theoretical 26.20 x 10% psi, experimental 
6 x< 105 psi and a present value of 24.35 x 10° psi. This 
data is summarized in Table V. x 


Reasonable correlation was also obtained in 


comparison with graphical results of Tsai [20] and Ekvall 16| . 


TABLE У. LONGITUDINAL STIFFNESS COMPARISON FOR THORNEL/EPOXY 


Longitudinal Stiffness, E. (psi х 107°) 


Analysis THORNEL 25* LTHORNEL 40** 
Blakslee [18] 12.78 26.20 
Experimental [18] 119397 24.60 
Present Analysis 12.44 24.35 


* THORNEL 25 Es/Ep 732 νεΞ.25 V 7.34 ve=.5 


** THORNEL 40 Ef/Ep=51 vg=.25 


45 





4. Composite Yield Strength 

Strength comparisons were more difficult to obtain. 
In some cases strength quantities described by identical 
terminology and applied to seemingly similar materials dif- 
fered by several hundred per cent. One reason for such 
discrepancies is the lack of consistent definition of 
strength terms. Valid comparisons only have meaning if 
the quantities investigated are precisely defined. 

Another important point is the effect of the matrix 
yield strength on composite yield strength. Matrix mate- 
rials of the same class, such as different brands of epoxy, 
normally have elastic constants that show little variation 
from one another. The difference in yield strength however 
can be considerable. The result of this difference is shown 
in Table VI, comparing a 40% filament model of Boron/2024 
Aluminum with a matrix yield of 1.30 x 10% psi and Boron/6061 
Aluminum with a matrix yield of 3.46 x 10* psi. 


TABLE VI. EFFECT OF MATRIX YIELD STRENGTH 
ON COMBOSTTE PROPERTIES 


νετ.» 
EC —  —""Boron/606l ^ Boron/2024  ————— — 
Aluminum . Aluminum $ Difference 
E; (psi)  20.53x10% 18.86x106 8. 
E, (psi)  31.07x10* 28.71x106 8. 
G (psi) 6.33x106 5.28х1 06 19. 
5 (051) 94.53х103 38.9x103 143. 
S, (psi) 28.50x103. 11.20x103 154. 
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As seen in Table VI the change in elastic properties 


is insignificant in comparison to the change in strength. 
The change that does occur in the elastic properties is due 
to minor differences in the elastic properties of the two 
matrix materials. The reason-for the very large increase in 


strength may be seen from Fig. (10). 





Figure (10). Strength Comparison of Different Alloy 
Matrix Materials. 


The slope of curve A represents the composite 
longitudinal stiffness, EM assumed to be identical for 
two composites, each having the same filament material 
but with different alloys of the same matrix material. 
As seen in Table VI the values of Er in such a case are 
quite similar. Curve B represents the stress-strain curve 
of one matrix material, curve C that of the other matrix 
material. They are assumed to have approximately the 
same elastic modulus (as, for instance, two different 
aluminum alloys) but with different yield stress. For 
the case of matrix B the yield stress will be reached at 


a composite strain Е The resultant composite stress 15 


ан 


сй... 
E 
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If the matrix material is that corresponding to curve 


C, matrix yield will occur at composite strain € resulting 


2! 
in composite stress S. . Thus a variation in matrix 
strength alone may markedly affect the composite strength. 

No comparison of strength terms should be attempted 
unless the complete properties of both constituents are 
specified. 

The example of Boron/2024 Aluminum of Table VI is 
directly comparable with the results of Adams 02. Рог 
identical material properties using a finite element 
analysis with constant strain triangles Adams reported a 
theoretical transverse yield strength of 11,200 psi, 
identical with the results of this analysis. 

Longitudinal strength may, be readily analyzed by 
strength-of-materials methods. Holister and Thomas [22 
developed a relationship for determining the ultimate 
tensile strength of a composite undergoing axial tension 


by considering the relation 
S “= бру + буш (III-13) 


where Ve and y are the filament and matrix volume fractions, 
respectively. 

For the case of axial load for continuous fibers, as in 
this study, the strain will be constant.  Holister and 
Thomas considered the particular case where the strain is 


such that the ultimate stress of the filament is attained 


while the matrix material is still in the elastic range. 
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This is equivalent to assuming a stress-strain relation 
corresponding to Fig. (11). Curve A represents the 
composite stress-strain relation and curve B that of the 


Matrix. 





Cult 


Figure (11). Stress-Strain Diagram for Failure Occurring 
in the Filament. 


If the matrix remains in the elastic region, the 


strain e€ for ultimate strength, S. must be less than 


ОЕ: Dt 
that necessary to cause initial yield in the matrix. Under 


these conditions 


n S ο Ξ On En (III-14) 
Equations (III-13) and (III-14) along with the relation 


+ у = l. (III-15) 
give 


Sult = Cult УҒЕр (1-Е /Е,) + £1] cÉm (III-16) 


Reidinger et al reported an equivalent form of this 


equation as 


SOV EE On (ITE 
ΠΠ: ° 





where Sm 15 defined as the stress in the matrix at the 


uit 
ultimate strain of the filament. This is shown in Fig. (11). 


Mié yield strength of the filament is Y, 

It must be repeated that these relations are for a 
composite that fails while the matrix is still elastic. It 
should be pointed out however that Holister pursues the 
matter further to consider the case of the matrix being in 
the plastic range. 

For the more common case of the matrix material 
yielding prior to the filament reaching yield, relations 
Similar to (III-13) and (III-14) may be utilized to find 
the initial yield strength. In the case of Fig. (12) the 
strain sufficient to cause initial yield for longitudinal 


load, SL is marked by ε 


D 





S 
L 
Yo 
ra 
Figure (12). Stress-Strain Diagram for Failure Occurring 


in the Matrix. 
For this case 


Y 
O O σ 
cu B — = Е CE 16 


Mere Yo is the axial yield of the matrix. Using Eqs. 
m 


(III-13) and (III-18) yields 


τ πο (III-19) 
m 
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The stress in the filament at the yield strain of the 


matrix 1s ge*, where 


s =“ v EE (III-20) 


The results for longitudinal strength calculated in 
this study differed by approximately 2% from those found 
by Eq. (III-20). 

5. Composite Poisson's Ratio 

а. Longitudinal Poisson's Ratio. The longitudinal 

Poisson's ratio, ντ, 38 calculated from the transversely 
isotropic stress-strain relations, Eq. (III-5), exhibits a 
linear behavior with respect to filament volume fraction. 
The results of Fig. (13) for E-glass/Epoxy correspond to 


those of a law of mixtures relation. 
Ур = УУ; + V_V (III-2]) 


This is in accord with the theoretical results of Rosen et 
al, 23] whitney and Riley, 17] and Tsai [20]. 

Specific experimental results were available from 
Blakslee et al [18] for THORNEL/Epoxy. This data also 
showed adherence to a rule of mixtures. 

b. Transverse Poisson's Ratio. The comparison of 


the transverse Poisson's ratio, with other theoretical 


Ут» 
studies raised several interesting points. Values found in 
the literature ranged from 0.7 29 ES considerably less 

кап 0.1. 18 AS the literature indicates a range ormene 


composite transverse Poisson's ratio that far exceeds that 
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of the Poisson's ratios of the constituents. The higher 
values exceed the theoretical maximum of 0.5 associated 
with isotropic materials. [24] 

In this investigation using the transversely isotropic 
stress-strain relations, Eqs. (III-5), results were obtained 
for ὯΝ With the addition of a small amount of filament 
material the composite transverse Poisson's ratio takes on 
a value greater than the Poisson's ratio of either constit- 
uent. For the materials considered the increase for a 30% 
filament volume was approximately 22% above that of the 
matrix. In all cases the value of Vin remained below 0.5. 

The transverse Poisson's ratio decreased with increasing 
filament volume. At high filament volume Vin values were 
less than the Poisson's ratio of either constituent. The 
maximum decrease noted in this study was 47% below the 
filament Poisson's ratio. 

To investigate the cause of the transverse composite 
Poisson's ratio having values outside the range of its 
constituents several additional computer analyses were 
made. The constituent and composite properties of these 
analyses are listed in Table VII. 

If the transverse Poisson's ratio was dependent only 
on the Poisson's ratios of the constituents and the filament 
volume, the results of analyses l, 2, and 3 below would be 
identical. As noted in the previous section, the longi- 
tudinal Poisson's ratio, u is dependent only on these 


quantities. This fact is borne out by the data of Table VII. 


53 





TABLE VII. TRANSVERSE POISSON'S RATIO STUDY 








Analysis E/E n Ve νε Vin ντ. Vin 
1 1 27 25 225 23 „250 
2 6 2 125 225 125 . 219 
3 120 el 25 25 225 109 
4 120 7 25 30 . 26 .134 


The results of this study indicate that the Poisson's 
ratio in the transverse direction is also greatly dependent 
on the ratio of constituent moduli of elasticity. The 
greater the ratio E C / Em! the greater the variance from the 
Poisson's ratio of the constituent materials. 

At small filament volume fractions the transverse 
Poisson's ratio exceeds that of either constituent. At 
large values of filament volume content Vip is considerably 
less than either constituent. At 50% filament volume the 
transverse Poisson's ratio is approximately equal to the 
average of the constituent Poisson's ratios, irrespective 
of the value of E /E This is shown in Fig. (14). 

Comparison of analyses 3 and 4 of Table VII shows 
that an increase in constituent Poisson's ratio does effect 
A but not as significantly as does the ratio E/E: A 205 
increase in v, resulted in approximately a 20$ increase in 


m 


Ут The effect of changing constituents Poisson's ratio 45 


limited however by the range of values for available 


materials. 
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migure (14). Effect of Constituent Stiffness on Transverse 
Poisson's Ratio. 


In the course of attempting to correlate the results 
of this study with other theoretical results several 
alternative analytical methods were considered. Some of 
these were found to give results for Poisson's ratio 
comparable to other extreme values found in the literature. 
These alternative methods were discarded as erroneous. To 
assist in evaluating the various predicted values of Ут 
two Of these methods will be discussed here. 

i) If it is assumed that the composite exhibits 
the same stiffness in all directions the stress-strain 


relations become 


Е = 5 (о - №0) 
x E x T y 
ze = 
Ey = (oy VmO,) 
ν 
a es 1 
€, ” ” — (o. + o.) + — O | (XII-22) 
E x y E Z 


Using the data of the various load cases these equations 
yield values for longitudinal Poisson's ratio, v y that 


follow a law of mixtures. As noted in the previous section 
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this concurs with most theoretical values. The transverse 
Poisson's ratio, however, takes on extremely small values, 
considerably less than .02 for 60% models. This method is 
not correct due to the inaccurate assumption regarding the 
stress-strain relations of the composite. 
ii) It was also attempted to evaluate the transverse 

ВЕ оп‘ = ratio by simulating a tension test. A finite 
element analysis was carried out on the model shown in Fig. 


(15a). The boundary conditions were 


u. = 0 on y= 0 


Y 
= 0 on x= 0 
= l on х= а 


This resulted in a displacement along the boundary as shown 


REIG. (15b). 


ү 





(a) Model 0 (b) Boundary Displacement 
Figure (15). Simulated Tension Test of a Composite. 
Consequently the Poisson's ratio, Vent 


boundary. In the problem considered the values ranged 


varys along the 


from approximately 0.15 to 0.7. The error in this method 


is that the rectangular boundaries of the basic block do 
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not remain rectangular, leading to violation of 
compatibility requirements. 

The above result does however suggest a physical 
understanding of the Poisson's ratio found in this 
investigation. Figure (16a) demonstrates how the results 
of this study require the composite to act under a tensile 
load. Figure (16b) shows how the material would react 
if compatible boundary conditions were not enforced. In 
Fig. (l6c) the stresses necessary to ensure the required 
compatible boundary conditions are shown. The stresses 
are such as to cause the displacements curve L to 
coincide with Lj: The integral of these stresses is zero, 
ї.е., S = 0. 


Y 
су? (on boundary) 0,=0 (on boundary) бу 70(оп boundary) 





(a) - 


Figure (16). Transverse Poisson's Ratio. 


Ec. Stress Concentration Factor (SCF) 

Before examining the results obtained for stress 
concentration factors some consideration must first be 
given to the value and interpretation of such results. 
The primary value of SCF is in its relationship to the 


strength of a material. For a homogeneous material 
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initial yield will occur at the point of highest micro- 
stress, that is the location of the SCF. For identical 
materials with different SCF, which may occur as a result 
of loading conditions, initial yield will occur at a lower 
stress level in the material with the larger stress con- 
Sentration factor. 

For composite materials SCF requires further interpre- 
tation. Due to material discontinuity, stress distributions 
will differ in the matrix and filament materials. For this 
reason SCF snould be calculated separately for filament and 
matrix regions of a composite. If the matrix is the larger 
no problem is posed, yield will occur in the matrix material. 
If the filament values are larger further investigation is 
necessary to ensure that initial yield is not occurring in 
the filament. If the ratio of filament to matrix stress 
concentration factor is less than the ratio of filament to 


matrix yield strength, yield will still occur in the matrix. 











That is 
Y 
SCF ¢ ος 
SCF ee] ушет in matrix 
m/ Yo | 
m 
SCF MO 
f f 
— = 1 simultaneous yield in 
SCF n ο matrix and filament 
SCF ¢ Yo 
zu 1 yield in filament (III-23) 
SGE ү 
m On 
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Caution should be exercised in interpreting stress 
concentration factors. Assuming that yield occurs in 
the matrix the macroscopic load required to cause yield 


may be related to the SCF as 


Sch = (III-24) 


where S. is the composite yield strength. This may represent 
a longitudinal, transverse, or longitudinal shear yield 
strength. For different materials with identical matrix 

E ud strength, Yo , therefore a smaller SCF will require 


In 


a larger macroscopic load, S to cause yield. Thus for 


ar 
such materials the SCF may be directly related to yield 
strength. An appealing aspect of this result is that for 
composites having identical matrix yield values, SCF offers 
an immediate evaluation of comparative yield strengths for 
the longitudinal, transverse, and longitudinal shear loading 
cases. 

For any composite with initial yield occurring in the 
matrix, regardless of the matrix yield strength, the compos- 
ite yield strength may be calculated for the three loading 


Cases as 


Sc = en (III-25) 


As noted above SCF also allows a rapid check on which 
constituent first reaches yield. 
For transverse and longitudinal shear loading the 


maximum SCF usually occurs in the matrix material. For 


59 





this reason any reference to SCF for these specific cases 
should be assumed to be a matrix property unless specifically 
stated to the contrary. No such generalization is warranted 
mee the longitudinal loading problem. 

Several authors have treated stress concentration factors 
extensively. Adams and Doner, [3] [6] using a finite difference 
technique, have provided a great deal of numerical data for 
both the transverse and longitudinal shear loading problems. 
For the transverse problem with a 55% Glass/Epoxy composite 
Adams and Doner calculated SCFX = 1.86. They also presented 
comparative data from solutions utilizing other analytical 
methods, ranging in value from 1.49 to 1.78. Analyses were 
not performed here for 55% volume composites, however а 47% 
model gave SCFX = 1.49, and a 60% model gave SCFX = 1.70. 


This data is included in Table VIII. 


TABLE VIII. SCF FOR TRANSVERSE NORMAL LOADING 





Method of Solution Stress Concentration Factor, SCFX 
Adams and Doner [6] 1.86 
Wilson [13] | a l 
Pickett [4] 1.49 
Foye [7] 1278 
Present analysis (473) 1.49 
Present analysis (60%) 1-70 
Present analysis (50% estimate) ο. 
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For the longitudinal shear problem Adams and Doner [3] 


gave results for SCFX of 40% and 55% models as 1.605 and 


1.787 for Glass/Epoxy and 1.649 and 1.854 for Boron/Epoxy. 


As shown in Table IX these values bracket the present 47% 


solutions of 1.66 for Glass/Epoxy and 1.69 for Boron/Epoxy. 


TABLE IX. 


Material 


Glass/Epoxy 
Glass/Epoxy 
Glass/Epoxy 
Boron/Epoxy 
Boron/Epoxy 


Boron/Epoxy 


40% 


4/5 


Filament Volume 


5СЕХ 


1.605 


1.787 


1.649 


1.854 


SCF FOR LONGITUDINAL SHEAR LOADING 


Source 


Adams and Doner 
Present Analysis 
Adams and Doner 
Adams and Doner 
Present Analysis 


Adams and Doner 
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IV. PARAMETRIC ANALYSIS OF COMPOSITES 


The purpose of this section is to analyze the effect 
that individual constituent material properties have on 
the mechanical properties of the composite. Properties 
considered include the modulus of elasticity, E, Poisson's 
ratio, v, and yield strength, E: of both matrix and 
filament materials. In this study it is assumed that both 
filament and matrix materials are isotropic. 

The goal is not to investigate every possible combination 
of material properties but rather to establish criteria and 
guidance for predicting composite properties associated with 
specific constituent properties. 

The macroscopic properties that are considered as per- 
formance indices include the initial yield strength, the 
specific strength, longitudinal and transverse stiffness, 
Shear modulus, stress concentration factors and specific 
moduli. 

To determine how constituent properties effect composite 
behavior, parameter variations were carried out on ficti- 
tious materials. Variations were performed by starting with 
Some actual properties, primarily those of a Boron/Aluminum 
or Graphite/Epoxy composite, and independently varying the 


particular property in question. 
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AN. EFFECT OF FILAMENT VOLUME 
To evaluate the effect of filament volume on composite 
properties several groups of 40%, 50%, and 60% filament 
volume models were investigated., These included 
1) Boron/NARMCO Epoxy 
2) High modulus graphite (MODMOR I)/4617 Epoxy 
3) High strength graphite (MODMOR II)/NARMCO Epoxy 
4) MODMOR 11/4617 Epoxy 
For MODMOR I/NARMCO Epoxy additional models of 30% and 70% 
were also utilized. The material properties for all mate- 
rials are listed in Table II. 
In a composite material subjected to a constant longitu- 
dinal strain the macroscopic stress is related to the micro- 


scopic axial stresses as 


Sc επ, - vor (IV-1) 


where Ve and Vin represent the volume fraction of filament 
and matrix materials. An assumption of this study is that 
the filament and matrix are perfectly bonded. For longitu- 


dinal loading the strain will be constant over any transverse 


cross section, thus; 


€ = Ef = € (IV-2) 
Within the elastic region the stress strain relations give 
Er, p UE C Y VER (IV-3) 


The longitudinal stiffness of a composite should be described 


then by a linear relation, the so called law of mixtures. 
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This line passes through the matrix modulus at 0% filament 


volume (E. = E, when Me = 0) and through the filament 


L 


modulus at 100% filament volume (E. = Е; when Ve = 1), as 
shown in Fig. (17) for MODMOR I/NARMCO Epoxy. 

For the case of transverse stiffness the relationship 
is not linear with respect to filament volume. For loading 
in the X or Y direction a constant strain applied at the 
boundary will not result in a uniform strain throughout the 
material. This is due to the change in material composition 
for cross sections perpendicular to the direction of loading. 
At zero filament volume the material has the properties of 
the matrix, that is the composite transverse modulus, Ent 
must equal the matrix modulus, Em. At 100% filament volume 
the material will have the properties of the filament and 


the composite transverse modulus, E will equal the filament 


πμ 
modulus, Er» The transverse stiffness, Em, versus filament 
volume, Ver curve must go through 
En = En at M 0 


and increase monotonically, terminating at 
En = Er at Vg = 1. 

All cases investigated followed this pattern. The results 
for MODMOR I/NARMCO Epoxy are plotted in Fig. (17). The 
increasing slope of the curve indicates that an increase in 
filament volume has a more significant effect on transverse 
stiffness at higher filament volume percentages. For example, 
for a MODMOR I/NARMCO Epoxy composite a 15% increase in 


filament volume, Ver from 35% to 50% results in a 42% increase 


in composite transverse stiffness, Ep. А similar 15% increase 
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Figure (17). Composite Elastic Properties vs. Filament Volume 
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in e from 60$ to 75$ results in a 75$ increase in composite 
transverse stiffness. For composite shear modulus, G, a 
change in cross sectional material properties indicates a 
behavior similar to that of the transverse stiffness. 

In longitudinal loading problem III, the stress con- 
centration factor SCFZ is calculated. The unit strain 
applied results in a uniform stress in the matrix material 
and a uniform stress of a different magnitude in the 
filament material. For the same constituent material and 
identical unit strain these stresses should not vary with 
changes in filament volume since they must satisfy the 
relations 


Оор ec (IV-4) 


lI 
tj 
m 


g (T75) 


πι πι 


The resultant surface traction, S5! the average value of 
the microstresses σε and a acting on the cross section, 
will change proportionally with the increase in filament 


volume, according to the equation 


S = VgO 


A AO (IV-6) 


£ 

Thus as (0,7) , the numerator of SCFZ, Eq. (III-1l), 
max 

remains constant, the denominator S, increases lineally 


with increasing filament volume, v and the stress con- 


Е 
centration factor itself decreases. 
For a homogeneous material the uniform stress throughout 


the entire block would equal the applied stress and SCFZ 
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would equal one. The curve for SCFZ versus filament 
volume should therefore pass through the points 
SOPZ (ve = 0) = 1, and SCFZ (У = 1) = 1. 

Consider the case of the SCFZ in the filament. The 
addition of a small amount of filament material to ἃ 05 
model would result in a negligible change in the surface 


fraction. That is, for unv+t strain 
є_ = 1.0 (IV-7) 


S; zes AE. E (IV-8) 


The maximum filament stress in this case would be 


σ Ε 

scrz = Í _ f πο 
σ E 
m m 


Figure (18) for Boron/6061 Aluminum shows this graphically. 
The points at 0% and 100% filament volume are established in 
accordance with the above discussion while the other points 
are results of the analysis. 

For the matrix SCFZ the maximum stress will remain 
constant while the average surface traction increases with 
increase of Ver as noted in Eq. (IV-6). For the matrix case 
therefore the curve will originate at SCFZ (νε = 0) = 1, 


and decrease to the point where, with little matrix material 


remaining, 


tx 


2 | 
D» (IV-11) 


SCFZ “ 





Recall from section II-B-6 that if the ratio of the 
filament SCF to the matrix SCF is less than the ratio of 
filament to matrix yield strengths, Y, s , yield will 
occur in the matrix. This will ον m each individual 
case which of the stress concentration factors discussed 
above is pertinent. 

SCFX and SCFXZ increase directly with filament volume, 
as shown in Fig. (18). These stress concentrations occur 


en the matrix at the interface. 


B. EFFECT OF FILAMENT YOUNG'S MODULUS 
TO evaluate the effect of variations in the modulus of 


elasticity of the filament, E on the composite mechanical 


к 
properties four computer analyses were made. Each analysis 
used a 50% filament volume model with material properties 


F Usted in Table X. 


TABLE X. FILAMENT YOUNG'S MODULUS VARIATION 


E ç O Em Om 


Analysis (рѕіх107%) о, (psixl0 °) (psixl0 °) о (psixl0”*) 


1 30 „2 4.8 10. e „3 
2 40 Д2 4.8 10. 3 $3 
3 ET . 4.8 10. 25. ‚3З 
4 70 τε 4.8 10. 2 58 
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Figure (18). Stress Concentration Factor vs. Filament 
Volune. 
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Figure (19) displays the transverse stiffness, longitu- 
dinal stiffness, and shear modulus of the composite as a 
function of fiber modulus. This figure shows that the most 
significant effect of increasing the filament modulus is 
the linear increase in the composite longitudinal stiffness, 


E This relates with matrix modulus variation and will be 


L' 
treated in greater detail in the following two sections. 
The relationship between the composite longitudinal 
stiffness, ы and the filament Young's modulus, Eg y is 
linear. In fact for the same reasons as discussed in 


section IV-A for filament volume variation, the relation- 


ship can be accurately described by a rule of mixtures as 


E, = νεξε + V E, (IV-12) 


The plot of the composite shear modulus variation of 
Fig. (19) appears to be linear, although it will be shown 
that this variation must be non-linear. Considerx that, for 
the combination of materials considered here, when the 


composite is homogeneous E, = 10 x 10° psi. Hence 


С = G, = G, Е ОИ 


The shear modulus curve of Fig. (19) must vary considerably 
from its apparent linear form to pass through this point, 
as shown in Fig. (25). 

As would be expected, considering the linear variation 
of longitudinal stiffness, the longitudinal yield strength 


also increased linearly with filament modulus increase. 
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Figure (19). 


Composite Elastic Properties vs. 
Filament Moáulus of Elasticity. 
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Variation in transverse and longitudinal shear strength 


was negligible for the range of E considered, Е; NUUS. X l0fpsi. 


C. EFFECT OF MATRIX YOUNG'S MODULUS 

To investigate the effect of variation of the matrix 
modulus of elasticity on the composite mechanical properties 
a series of four analyses were made. All analyses were made 
with a 50% filament volume model with constituent properties 


as shown in Table XI. 


TABLE XI. MATRIX YOUNG'S MODULUS VARIATION 


Er ο J Em MOL 
Analysis  (psix10 5) v,  (psixl0 ?)  (psixlO °) v,  (psixl0”*) 
1 60. „2 4.8 10. . . 346 
2 60. ο 4.8 T 3 .346 
3 60. 2 4.8 Ir „3 ‚346 
4 60. 2 4.8 2.5 5 .346 


As established in Eq. (IV-3) the composite longitudinal 
stiffness varys lineally as 


Bro = УЕ. + у E 
L ΓΕ m m 


Thus the variation of Er with matrix modulus is also linear 
as shown in Fig. (20). From this figure 1t is seen that the 
variation of E, has a greater effect on the transverse 
stiffness and shear modulus than on longitudinal stiffness. 
An increase of 400% in the matrix modulus of elasticity 


results in percentage increases of 
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Figure (20). Composite Elastic Properties vs. 
Metrix Modulus of Elasticity. 
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iii) G = 222% 

This also resulted in an increase in transverse and longitu- 
dinal shear yield strength. The complementary effect of 
this result with that of the filament modulus variation 
will be discussed in the following section. 

An interesting aspect of the variation of the matrix 
Young's modulus is that, despite the increase in longitudinal 
stiffness, the longitudinal yield strength, Si; decreases. 


The reason for this may be seen in Fig. (21). 


S — ome на алә «ο «| απ == 
L, | 
| 
S ———— | 
L | 
t | | 
I | 
ον ο τς = 
m | 
| | 
| 1 
E 
ee | 
Eure (21). Effect of Matrix Stiffness on Longitudinal 


Yield Strengtn. 


For two matrix materials with similar initial yield points 
one with the lower Young's modulus will require a larger 

Strain to reach yield. Due to the relatively small effect 
of matrix stiffness on composite longitudinal stiffness, Er, 
assume that E, is the same for both cases. The larger strain, 


£5; will result in a larger composite yield strength, S; . 
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BESZFERECT OF CONSTITUENT PROPERTY RATIO 

The results noted for filament and matrix stiffness 
variation may be conveniently correlated through use of 
non-dimensional parameters. For example, composite longi- 
tudinal stiffness, E may be normalized with respect to 
filament stiffness, i.e., Ετ/Ες, and plotted against the 
constituent modulus ratio En Eg - Figure (22) is such a 


representation. Longitudinal stiffness may be expressed as 


Dividing through by Ef and rearranging terms yields 
E, /Eg * Vg (1 - En/Eg) + ЕЕ; (IV-14) 


Equation (IV-14) describes the configuration of the non- 
dimensional longitudinal stiffness curve. The following 
points are significant: 


i) when E. = Ec, i.e., a homogeneous material, 


E, /Ef = En Ee = 1. (IV-15) 
ii) when pm = 0, 1.e., zero matrix stiffness, 
E, /Ef = УЕ (1) = e (IV-15) 


Equations (IV-15) and (IV-16) were used to determine the 
extreme points of Fig. (22). The points noted on this 
figure are the results of this investigation for a filament 
volume fraction, Ve = .5. For comparison some other curves 


were drawn using the endpoints determined by Eqs. (IV-15) and 
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(IV-16). The non-dimensional longitudinal stiffness 
curve for any filament volume must fall in the triangular 
Region of Fig. (22) formed by the vs = 0 curve and the 
ИЕ = 1.0 curve. 

The comparative effect of filament Young's modulus 
and matrix Young's modulus een on composite longi- 
tudinal stiffness may be examined by consideration of Eg. 
(IV-13). Taking partial derivatives of this equation with 


respect to the individual independent quantities yields 


L 
== y (IV-17) 
JE i 
m 
and 
JE 
L 
JE; 


Thus the effect of constituent stiffness variation on longi- 
tudinal stiffness is directly related to the volume fraction 
of the constituent. For the 50% filament volume case this 
is seen to be true as Figs. (19) and (20) show that filament 
and matrix stiffness variation, respectively, result in 
identical slopes of 0.5 when plotted versus longitudinal 
stiffness. 

These comments concerning comparative effects of matrix 
and filament stiffness variation refer to the change per unit 
stiffness of constituent. The actual range of longitudinal 
stiffness due to matrix and filament stiffness variation is 
considerably different however. Consider Fig. (23), a non- 
dimensional longitudinal stiffness curve for some arbitrary 


filament volume fraction, νε. 
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Ereoer of Constituent Stiffness on 
Longitudinal Stiffness. 
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Figure (23). Non-Dimensional Longitudinal Stiffness. 


For a given constituent moduli ratio, say point A, a 


specific ratio, E-/Egr results at point A'. At this point 


If the constituent ratio is increased to point B by increas- 
ing the matrix modulus while holding the filament modulus 


constant then 


tr 

1 

tJ 

ty 
Γη 


The ratio of longitudinal stiffness for these two hypothetical 


composites is 


t 


E TAE = АТ TEEN p. 
A A A 


Thus variation in the matrix modulus results in a change in 


composite longitudinal stiffness in direct proportion to the 
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change in non-dimensional stiffness, E, /Eg- In general 
the maximum range of non-dimensional stiffness, E/Egr is 
from the value of filament volume fraction, Ver tOMPSOS 
The maximum increase in longitudinal stiffness Бас due 


solely to a variation of matrix modulus is 


E 
Le 1.0 


ΠῚ -π 
En t 
max 


For the 50% filament volume case of Fig. (22) this maximum 


Encrease is a factor of 


B 1.0 


LE SE 80 
L E 
max 


If the constituent ratio of Fig. (23) is increased to 
point B by decreasing the filament modulus while holding 


the matrix modulus constant, then 


Е а (1ν-20) 
Thus the increase in longitudinal stiffness with variation 


of the filament modulus is dependent on the non-dimensional 


stiffness, E/Egr аша пе ratio of filament moduli, ὃς Ve: 
A B 
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For actual filament materials the range for the filament 
modulus is approximately 10 x 10% to 60 x 10% psi. Consider 
a constant matrix material, for example aluminum. For the 
50% filament volume case of Fig. (22) the range of E, /Ef 
corresponding to this range of filament modulus is 50 to 
1.0. For this 50% case the ue dion in longitudinal stiff- 


ness due to variation in the-filament modulus alone is 


Er, 


> 


_[.50) [so 


Er, B1 919 
B max 


= 3,0 


Variation of the filament modulus is thus seen to result in 
a much greater range of longitudinal stiffness than does 
the variation of the matrix modulus of elasticity because 
of the greater range of Ef as compared to Em" 

Figure (24) is a non-dimensional representation of the 
effect of constituent stiffness, En/Eg y on a normalized 
composite transverse stiffness, E Eg. Note that the curve 
approaches Em = 0 as En Ef approaches zero. This indicates 
that when the matrix material is extremely soft in comparison 
to the filament that the macroscopic transverse stress 
required to produce a unit displacement of the block boundary, 
as in load case I, 15 virtually zero. In the limit, when 
En = 0, then E, = 0. Thus it is theoretically possible to 
produce a composite with a value of transverse stiffness 
approaching zero. 


Comparison of Figs. (19) and (20) shows that the slope 


of the En vs. matrix stiffness curve is greater than that 
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Figure (24). Effect of Constituent Stiffness on 
Transverse Stiffness. 
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} 


of the En vs. filament stiffness curve. Therefore on a 

unit stiffness basis, variation of matrix Young's modulus 

is more effective than variation of filament Young's modulus 
on the composite transverse stiffness. 

By considering a non-dimensional transverse stiffness 
νε and following the same technique as for the non- 
dimensional longitudinal stiffness the range of the trans- 
verse stiffness due to matrix Young's modulus variation, 
analogous to Eq. (IV-19), is 

Em (Ep/E£) 


—B n (IV-21) 


E 
T (Ep/ Ec) 
max max 
The range of E Eg as is seen in Fig. (24), is zero to one. 
Thus the theoretical maximum increase of transverse stiff- 


ness due to matrix Young's increase is 


B de 0 


ЕЕЕ ER a = 00 


E = = 


max 


The maximum variation in transverse modulus due to variation 


of the filament stiffness alone is 


Ep. (Б/Е) Eg 
= = sacs alls (IV-22) 
E (ZS Er 
T quaest 
B| max B В шах 


For common filament materials the range of filament modulus 
is approximately 10 x 10° to 60 x 10% psi. Consider a 


constant matrix material, for example aluminum. For the 50% 
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filament volume case of Fig. (24) the range of E Pe 
corresponding to this range of filament modulus is .349 
to 1.0. Thus the maximum variation in transverse stiff- 


ness due to variation of the filament modulus alone is 


Е 
Tal  [.349| 60. 


Í 3 — 
x. 


Е ДЕт` ο. τας. 
B/ max 


= 2.094 


The range of transverse stiffness due to matrix Young's 
modulus variation is greater therefore than that due to 
filament Young's modulus variation. 

Figure (25) is a plot of the normalized composite shear 
modulus, G/G., versus the ratio of constituent shear moduli, 
G Gç. This curve substantiates earlier contentions that 
shear modulus variation with respect to constituent moduli 
is non-linear. Again this curve approaches zero as the 
constituent moduli ratio approaches zero. 

Again comparing Figs. (19) and (20) it is seen that 
the slope of the composite shear modulus versus matrix 
stiffness is greater than the slope of the composite shear 
modulus versus filament stiffness curve. This indicates 
that, per unit of constituent stiffness, the matrix material 
has a greater effect on composite shear modulus. 

By considering a non-dimensional shear modulus curve 


and following the same technique as for the non-dimensional 


longitudinal stiffness, the range of shear modulus variation 


due to matrix shear modulus variation, analogous to Eq. (IV-19), 


is 
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(С/С) 
Gp f = 
= = (6787 (IV-23) 


А | max Е A | max 


The range of G/G., as seen in Fig. (25), is from zero to 


1.0. Thus the theoretical maximum value of Eq. (IV-23) is 


р К В TO 
0 


A max 


If the ratio δε 15 increased by changing the filament 


modulus the maximum range attainable is 


“a = E “E, 
SE max (G/G£) Gf 
B В | пах 


(IV-24) 


For common filament materials the range of shear modulus is 
approximately 4.0 x 10° to 25.0 x 10% psi. For a constant 
matrix material, for example aluminum, the range of G/Ge ОЕ 
Fig. (25) corresponding to this range of filament shear 

modulus is .31 to 1.0. These values are for a 50$ filament 


volume case. For this case the Eq. (IV-24) becomes 


G 

_А s «31 25 1.94 

Gp IO 4 ° 
max 


Variation of the matrix shear modulus thus results in 
a greater range of composite shear modulus than does 
variation of filament shear modulus. 

The effect of Poisson's ratio will be discussed in the 


next section. 
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Figure (25). Effect of Constituent Shear Modulus on 
Composite Shear Modulus. 
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The primary advantage of the non-dimensional represen- 
tations is that valid comparisons may be made between 
composites without necessarily holding the Young's modulus 
of one constituent constant. 

It should be noted that the non-dimensional plots 
presented here are for a 50% filament volume model. Dif- 
ferent A o e models will result in different 


Curves. 


e EFFECT OF POISSON'S RATIO 

Two series of computer analyses were made for evaluation 
of the effects on composite behavior due to change of Poisson's 
ratio of the filament, v;. Four analyses were made using a 
Fictitious material approximating the properties of Boron/ 
Aluminum while varying the filament Poisson's ratio from .15 
to .30 in increments of .05. 

Another series of analyses was conducted with MODMOR I 
graphite and epoxy. This was done because Poisson's ratio 
for graphite filaments is not accurately specified and it 
was desired to determine whether or not this would signif- 
icantly effect the results of this investigation. 

In both cases no significant effect was noted on any 
composite mechanical property. All observed variations 
were less than 13. 

The effect of the variation of matrix Poisson's ratio 
was evaluated through a series of three analyses. Gonstituent 
properties similar to Boron/6061 Aluminum were used while 


varying the matrix Poisson's ratio from to .30 to .35. Again 
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no significant effect on composite mechanical properties 
was noted. It is expected however that changes in Poisson's 
ratio will produce more significant changes for post yield 


behavior. 


КОО EFFECT OF MATRIX YIELD STRENGTH 


ШОШ сс солгу the effect of the matrix yield strength, 


ES, on composite mechanical properties, a series of four 


m 
computer analyses were made. All used a 50% Boron filament 


volume model with constituent properties as listed in Table 


XII. 


DABLE XII. MATRIX YIELD STRENGTH VARIATION 


Е Yo. E. Yo 
Analysis  (psixlO 5) Ve (pš ix]0 °) (psixl107*) ee (psi 10m) 
1 60. 2 4.8 10; τη . 346 
2 GNE -ο 4.8 πο ο. 

3 60. 2 4.8 ΠΠ 3 «150 
4 60. 2 4.8 ΠΟ 23 . 030 


The yield point of a material has no effect on the 
elastic properties of that material. It merely defines the 
limit of elastic behavior. Thus the variation of matrix 
yield strength had no effect on the composite stiffness or 
Shear modulus. This is as expected considering the compar- 
ison in section III-B-4 of two aluminum alloy matrix 


materials. 
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Figure (26) shows hypothetical stress-strain curves for 
the composite (curve A) and the matrix (curve B). For a 
given matrix yield strength, say ү , a Specific composite 
macroscopic yield strength, Sor is defined. If, as stated, 
the yield point of a material does not effect the elastic 
properties, then the stress-strain curves remain the same, 
independent acne yield strength variation. For a given 
increase in matrix yield strength, say a 50$ increase to 
M , plane geometry shows that an identical increase, on 

m 


a percentage basis, must occur in the composite yield 


strength. Figure (27) confirms this hypothesis. 





Figure (26). Effect of Matrix Yield Strength Variation. 


It should be noted that the above discussion assumes 
that yield occurs in the matrix, as is true for most, but 
not all, composites. 

The implication of this phenomenon is that there is no 


variation in the shape of the initial yield surface due 
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solely to variation in the matrix yield strength. That is, 
the initial yield surface changes in size only. This can 


be seen in Figs. (C-1) and (C-2) of Appendix C. 


G. EFFECT OF FILAMENT CROSS SECTION 

The majority of the computations included in this investi- 
gation are based on longitudinal filaments with a circular 
cross section. This was done because it is assumed that this 
particular cross section would be the easiest to produce. The 
question arises as to the effect change of shape of filament 
cross section would have on composite properties. Such 
changes might be the result of fabrication. 

To provide an answer to this question two elliptical 
cross sectional filaments were investigated. These were 
identical ellipses with major/minor axes ratio of 1.25 to 
1.0. One model had the major axis oriented along the x-axis 


and one along the y-axis, as seen in Fig. (28). 


x axis Major Axis y axis 


πε (25). Bide tical Cross Sections. 


The rule of mixtures and the comments of section IV-A 
leading to its use are independent of cross section shape. 
Elliptical and circular cross section filaments of the same 


volume should therefore produce identical results for the 
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longitudinal composite stiffness, E, > In the three cases 
investigated, 50$ models of Boron/6061 Aluminum, Boron/NARMCO 
Epoxy, and S-glass/NARMCO Epoxy, this was found to be true. 
The variation in E, was less than 1% in all cases, as may be 
seen from Table XIII. 


- TABLE XIII. EFFECT OF -FILAMENT CROSS SECTION 
ON LONGITUDINAL STRENGTH AND STIFFNESS 


Circular Elliptical 
. "6 » u » = 6 ۰ “h 
Е; (р51х10 ) Жел ) E, (psixl0 ) S, (psix10 ) 


Boron/6061 Aluminum 37.08 11.4 37225 102: 
Boron/NARMCO Epoxy 287229 1679 2/9857 14.5 
S-glass/NARMCO Epoxy 6.2.05 3. 47 6.105 3.04 


The transverse composite stiffness reacted as expected. 
In all cases En increased in the direction of the major el- 
liptical axis and decreased in the direction of the minor 
axis as noted in Table XIV. This seems reasonable since more 
of the stiff filament material is oriented in the direction 
of the major axis. Оп а percentage basis the increase in 
transverse stiffness in the direction of the major axis was 
approximately twice that of the decrease in the direction of 
the minor axis. 

An interesting result of the cross section variation 
that is evident in Table XIV is that the transverse strength 
changed in a manner ee to that o£ the transverse stiff- 


ness. That is, the transverse strength increased in the 


direction of the minor axis and decreased in the direction 
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of the major axis. This result may be understood by 
considering the micromechanics of the problems in question. 
For ease of discussion call the transverse load case with 
the major axis oriented along the x-axis load case IA, and 
that with the major axis along the y-axis as load case IB. 
These load cases are shown in Fig. (29). Using this 


notation, the analysis results require, for transverse 


normal loading, at yield 
ο (IV-25) 


The definition of macrostresses requires that the integral 


of the microstresses along the boundary be related as 


E [οι (1ν- 26) 


TOBLE 215. ΕΟΤ OF FILAMENT CROSS SECTION 
ON TRANSVERSE STRENGTH AND STIFFNESS 


circular Fili I ico] 
lajor Axis Minor Axis 


En Sm En Sm 











Boron/6061 Aluminum 
Boron/NARMCO Epoxy 
S-glass/NARMCO Epoxy 


F (psix10°) 


T ST (р51х10*) 


Figure (29) shows the stress and strain distribution on 


the x-boundary (x = a) for a unit displacement for load 
cases IA and IB. In both cases yield occurs in the matrix 


at a point near the interface, along the x-axis. Both cases 
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also show little variation in stress or strain along the 
x-axis between the interface and the boundary. For example, 
for load case A, i.e., major axis along the x-axis, problem 
I for a 50% S-glass/Epoxy model gives a stress, o,, at the 
interface of .664 x 10° psi. At the boundary (x = a) o, 
equals .688 x 10% psi. The same material for load case B 

gave an ai cr Eê value of o, as .366 x 10° psi and a 
boundary value of .389 x 10° psi. Thus the strain to cause 
yield will be approximately the same at the yield point, i.e., 
the interface, and at the boundary (x = a, y = 0). If the 
Strain and stress curves of Fig. (29) are scaled to give 
identical strain values at the boundary (x = a,y = 0) the 
more uniform distribution of case IB will result ina larger 
value of the integral of the stress. This confirms Eqs. 
(IV-25) and(IV-26). 

It seems therefore that the stiffness of the composite 
depends on the degree of filament material alignment in a 
particular direction while the strength depends on the 
maximum filament cross section in a direction perpendicular 
to that of the load. The inference of this is that the 
optimum cross-section shape would approximate that of a 
plane whose orientation depends on the desired properties. 
This is shown in Fig. (30). The filament shape can not be 
a rectangle as the square ends would induce prohibitive 


stress concentrations. 
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Y os 
x 3 
For Maximum E, For Maximum Sy 
Figure (30). Approximate Optimum Cross Section for 


Particular Properties. 


In the case of circular cross section filaments the 
transverse properties are the same in the y-direction as 
in the x-direction. Thus transverse yield strength in 
the y-direction would be numerically equivalent to that in 
the x-direction. This is shown in Fig. (3la). Again assume 
the convenient designation of load case A as that of the 
elliptical cross section with the major axis along the 
x-axis, and load case B with the major axis along the y- 
axis. The results of these two load cases may be interpreted 
as transverse strength in different directions of the same 
model as shown in Fig. (31b). 

On a percentage basis the change in transverse strength 
due to cross section variation was considerable, ranging from 
153 for Boron/Aluminum to 23$ for Boron/NARMCO Epoxy. O£ 
more significance however is the fact that the vector sum 
of S. and Ὃν is independent of the filament cross section. 


This indicates, in vector notation, for Fig. (31) 


SS | (IV 27) 
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Or, alternatively, 


DNE εις 5] (1ν-- 28) 





Y Circular Filament 
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Elliptical Filament, Identical v. 
(b) 
agure (31). Transverse Strencth Relation to Filament 


Croesi Section: 


This indicates that any increase in strength in one direction 
is directly offset by a decrease of the strength in the per- 
pendicular direction. 

In each case there was a noticeable increase in the stress 
concentration factors SCFX and SCFXZ, although SCFZ remained 
relatively constant. The net result was a reduction in the 


longitudinal yield strength of the elliptical cross section 
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compared to the circular cross section. The decrease was 
approximately 12% in each case. The implication is that a 
circular cross section is the most effective, with respect 
to a given filament volume, in terms of longitudinal strength. 
Adams and Doner [6] meperted tNat stress concentrationm factors 
are directly related to the maximum radius of curvature of 
the ment. For a given filament volume a circular cross 
section results in the smallest possible maximum radius of 
curvature. This also implies that the circular cross section 
is the most efficient. 

Another significant result evident in the analysis of 
the elliptical filament is the effect of filament shape on 
the composite shear modulus, С. Normalized with respect to 
the circular composite shear modulus the results obtained 
are listed in Table XV. 


TABLE XV. EFFECT OF CROSS SECTION VARIATION 
ON COMPOSITE SHEAR MODULUS 


v= .5 
Ge11/G8c f 


Composite Major Axis Orientation 


Boron/6061 Aluminum 5 $4555 
Boron/NARMCO Epoxy d] . 326 
S-glass/NARMCO Epoxy T15 .398 


The shape of the filament cross section could be used 


to good advantage by a designer with detailed knowledge of 


97 





the expected loading of the structure under design. Such 
large variations in some composite properties due to 
relatively small departures from a circular cross section 
also make it extremely dangerous to load a "circular" 
filament near its expected transverse or shear strength. 
Small manufacturing errors can lead to large strength 


discrepancies.-- == 
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V. COMPARISON OF MATERIAL COMBINATIONS 


The preceeding sections have investigated the effect 
of various constituent properties on the properties of 
the composite. The goal was to be able to predict the 
Eemparative properties of various actual material combina- 
pons. Іп ТЕ section several functional material com- 
binations are analyzed. Theoretical mechanical properties 
are compared and the results are correlated with the guide 
lines developed in the preceeding sections. 
Forty, fifty, and sixty percent filament volume models 
are compared for the following combinations of materials, 
i) MODMOR II/NARMCO Epoxy 
ii) MODMOR II/4617 Epoxy 
111) MODMOR I/4617 Epoxy 
iv) Boron/6061 Aluminum 
For MODMOR I/NARMCO Epoxy and E-glass/Epoxy, 30$ and 70$ 
filament volume models were also considered. Single models 
of S-glass/NARMCO Epoxy, THORNEL 25/4617 Epoxy, THORNEL 
40/4617 Epoxy, and Boron/2024 Aluminum were also investigated. 
Primary measures of structural efficiency of a material 
are the specific strength and specific moduli of the material. 
These are defined in section I-B, however their importance 


necessitate further explanation. 


In general, design of a particular structure requires 
certain material properties of strength and/or stiffness. 


in modern technology the weight of a structure has gained 
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equal importance with mechanical properties. It is 
significant to note that all common structural materials, 
when considered on a stiffness to weight basis, are remark- 
ably similar. Wood, steel, aluminum, and titanium fall 
hin the range of 92.0 x 10 f psi/(lb/in?) to 

105.5 x 10° psi/(lb/in’) for specific Young's modulus. 
u for specific ultimate strength these common 
Materials range from 4.60 x 10° to 8.33 x 10° psi/(lb/in’). 
This data is summarized in Table XVI. Thus for a composite 
to be valuable from a structural standpoint it must exceed 
these values for specific mechanical properties. 


ΠΕΙ XVI. SPECIFIC ULTIMATE STRENGTH AND SPECIFIC 
YOUNG'S MODULUS FOR COMMON STRUCTURAL MATERIALS 


Specific Ultimate Strength Specific Modulus 


EN erial pei (Ipin )xlo ° psu Ab din sto"? 
Steel 5.16 105.5 
Aluminum 6.64 104.2 
Titanium S53 9352 
Magnesium S 101.6 
Wood 4.60 92.0 


Flgure (32) shows the transverse stiffness, Ent and 
the specific transverse stiffness, En/ pr where p is the 
composite specific weight. These results are for the 50% 
filament volume models considered. In the preceeding 
sections transverse stiffness was found to depend primarily 


on the matrix modulus of elasticity. The results of Fig. 
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(328) show that the transverse stiffness of materials with 
similar matrix properties showed little variation. For 
example, those composites with a matrix of NARMCO epoxy 
had comparable values of transverse stiffness regardless 
of filament properties. The aluminum matrix was so much 
stiffer than any other matrix material studied that the 
performance ren ken far exceeds that of the 
other composites. 

Figure (32b) shows the specific transverse stiffness. 
Note that the smallest value of composite specific transverse 
stiffness is approximately 80% greater than that normally 
associated with structural materials. 

Figure (33) shows the composite transverse strength, Sp’ 
and specific transverse strength, S/P- The constituent 
property that showed the most influence on transverse 
strength was the matrix modulus of elasticity. Figure (33a) 
reinforces this judgement by showing negligible variation 
in transverse strength among those composites utilizing 
identical matrix materials. Progressively stiffer matrix 
materials produced greater transverse yield strength regard- 
less of the filament material. 

Figure (33b) shows the specific transverse yield 
strength for the composites considered. Note that not 
all of the material combinations exceed the stated value 
of 8.33 x 10°in. for common structural materials. One 


reason for this is that the quoted values for common 


structural materials are ultimate strengths rather than 
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Figure (32). Transverse Stiffness Comparison. 
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Figure (33). Transverse Strength Comparison. 
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yield strengths. Another reason is that the composite 
lamina are designed such that their primary strength is 

in the longitudinal direction whereas most standard struc- 
tural materials are considered isotropic. 

Figure (34) is a plot of longitudinal stiffness and 
specific longitudinal stiffness. It has been shown in 
this НЫЕ composite longitudinal stiffness can be 
accurately predicted by a law of mixtures. For the 50% 
case of Fig. (34a) the composites should be in order of 
increasing sum of matrix and filament moduli of elasticity. 
This is seen to be true. 

Figure (34b) shows the specific longitudinal stiffness. 
Note that the graphite/epoxy combinations are considerably 
more efficient than the Boron/Aluminum combinations despite 
the large stiffness advantage of aluminum over epoxy. 

Figure (35) shows the longitudinal yield cur аша 
specific longitudinal yield strength of the various 
composites. 

The relative longitudinal yield strength of a composite 
is somewhat more difficult to predict than are the other 
properties. Increasing the filament modulus or the matrix 
yield strength has been shown to increase the composite 
longitudinal yield strength. Increasing the matrix modulus 
of elasticity decreases longitudinal yield strength. _ There 
is no set method however for predicting the combined effect 


ӨТ sueh changes. 
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Figure (34). Longitudinal Stiffness Comparison. 
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Pigure (35). Longitudinal Strength Comparison. 
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The MODMOR II (graphite)/Epoxy models indicate a 
case where the matrix yield strength was of greater 
significance than the matrix modulus increase. The 4617 
Epoxy had a greater modulus, which should decrease the 
yield strength, and a greater matrix yield strength which 
should increase the composite longitudinal strength. The 
net а marked increase in composite yield 
strength. Considering the variations individually, as 
in Fig. (36), this indicates that the increase in matrix 


yield strength, say from Y, to Y, ^, increased the composite 


m m 
yield strength from s. to Sr. This is shown in Fig. (36). 


The decrease in matrix modulus results in a composite stress- 
strain curve with a smaller slope. For the same matrix 
yield strength, Y. °, the composite yield is ST. The results 


m 
of the analysis shows that 


Sr < ex < S (V-1) 





Ë 


€ 


О 





(b) 


Figure (36). Combined Effect of "Matrix Modulus and 
Matrix Yield (MODMOR II/Epoxy) 
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The two boron analyses represent the opposite case. 

The change in matrix modulus prevails over the change in 
matrix yield strength as the epoxy composite has a yield 
strength considerably larger than that of the Boron/Aluminum 
composite. 

It should be noted that the longitudinal yield strength 
of the К карт се composite MODMOR 1/4617 Epoxy is considerably 
less than that of the MODMOR 11/4617 Epoxy. This is signif- 
icant in that the Young's modulus of MODMOR I is 50$ greater 
than that of MODMOR II. This apparently contradicts the 
guidelines presented in this study. The reason for this is 
that the conclusions of this investigation are based on the 
assumption that initial yield occurs in the matrix rather 
than the filament. Due to the relatively low yield strength 
of MODMOR I initial yield actually occurs in the filament. 
This is predicted by considering the stress concentration 
factors as outlined in section III-B-6.  Explicitly, for 
MODMOR I/4617 Epoxy, 


ΕΙ ποιο «δια νο αμ 2-5 x 10% psi 


Matrix yield strength -08 x10 psi 
Filament SCFZ 2.106 
Matrix SCFZ .0524 


SCF2g 2.106 








SCFZ ^ .0524 5 40.2 (У-2) 
Уо 
Bus 
To = 08 = 31.25 (V-3) 
m 
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Which meets the criteria of section III-B-6 for initial 
yield occurring in the filament. 

Figure (35b) represents the specific longitudinal yield 
strength. . The unshaded area.of the bar chart represents 
the specific yield strength from the yield calculations of 
this analysis. The shaded region of the figure represents 
an estimate of specific ultimate strength. The ultimate 
strength was estimated by use of the equation 


S = V Y +v. o (V-5) 


As discussed in section III-B-4 E is the stress in the 
ult 


matrix at the ultimate strain of the filament. Again it 
should be noted that this relation assumes that the matrix 
remains elastic up to the failure of the filament. Generally 
this is not true. Using the method of Lin, Salinas, and 

Ito [ο] the theoretical ultimate tensile strength of a 50% 
Boron/Aluminum composite was found to be 243,000 psi. This 
compares with an estimate using Eq. (V-5) of 305,000 psi, a 
difference of greater than 20%. While the discrepancy is 
rather large the results are considered to be worthy of 
inclusion as a relative measure of the range between yield 


and failure. 


For MODMOR I, MODMOR II/4617 Epoxy and THORNEL 25 (all 


graphite filaments) yield actually occurred in the filament. 
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This is the reason that the specific ultimate and specific 
yield strength are virtually coincident. 

Also of importance in a manufacturing process is the 
cost of the materials. Table XVIII shows strength per 
unit cost data for selected materials. Cost data was 
obtained from Refs. [25] and [26]. The extremely high 
strength to cost ratio for glass filament serves to explain 
why glass composites are in such wide use despite their 
relatively low transverse yield strength. The present 
interest in graphite filaments can certainly be justified 


by the strength/cost increase anticipated in the near future. 


TABLE XVII. COMPARATIVE STRENGTH PER UNIT COST 
Epoxy Matrix Ve = 0.5 
Filament Strength/Unit Cost (in-lb/dollarx10 5) 
Boron . 98 
Graphite ZERO 
Graphite (1974) 4.02 
Glass 42.48 


Numerous other factors also enter into a material 
selection process. Ease of handling and fabrication, 
Special processing requirements, and specific desirable 
properties such as thermal expansion may be such factors. 
Consideration of these factors is beyond the scope of this 


investigation. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

In this investigation a finite element analysis of the 
microscopic and macroscopic stresses of a composite material 
has been carried out. The stresses were used to calculate 
the initial yield surface for unidirectional composite 
materials. The results of these analyses were then used 
to determine the mechanical properties of a laminate of 
composite material. Parametric studies of the effect of 
various constituent properties on the composite mechanical 
properties were then carried out. Table XVIII summarizes 
the results of these parameter variations. 

These results are based on the assumption that initial 
yield takes place in the matrix material. The effects noted 
in this table are those associated with a change of the 
specified constituent property only, that is, with all 
other constituent properties held constant. Caution must 
be exercised in attempting to predict the effects of changes 
ln more than one constituent property. 

Other results included: 

1) Composite Transverse Poisson's Ratio, Vint This 
value exceeded that of either constituent. By adjusting 
the ratio of filament Young's modulus to matrix Young's 
modulus, E / Eg the designer may achieve extreme values of 


transverse Poisson's ratio. 
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ii) Longitudinal Initial Yield Strength, Spy. For 
initial yield occurring in the matrix material the longi- 
tudinal yield strength may be predicted by 


E = У СЕ + p 


where е is the yield strength of the matrix material, 
of is ee in the filament at the yield strain of 
the matrix, and νε and v, are filament and matrix volume 
fractions, respectively. 

iii) Initial Yield Location. The computer program 
utilized in this investigation gives results which include 
the location of the yield point for any macroscopic stress 
state. For load case III it may be determined whether the 
composite fails in the matrix or the filament by r na s 
ing the stress concentration factor in the longitudinal 
direction, SCFZ. If the ratio of SCFZ in the filament to 
that in the matrix is ied than the ratio of filament yield 
strength to matrix yield strength then yield will occur in 
bue matrix. That is, if 

SCFZ¿/SCFZ, < Yo /Y, 
Е m 
yield will occur in the matrix. This determines the appli- 
cability of the longitudinal yield equation noted above and 
does not require the application of any particular yield 
criteria. 
iv) Filament Cross Section Effects. Variations from a 
clrcular filament еве section have significant influence 


on composite properties. Shear modulus and transverse 


Ti 





properties change considerably. Due to the distribution 
of microscopic stresses on the block boundaries an el- 
liptical cross section has greater transverse strength in 


mae direction of the minor axis. 


B. RECOMMENDATIONS 

i) The logical extension aE this investigation is 
ee ideration of inelastic loading. The analysis used 
in this study is a special case of the elastic-plastic 
analysis of Lin, Salinas, and Ito. [ο] 

ii) Further effort seems warranted in the evaluation 
of filament cross section effects. Through use of cross 
sectional shapes and filament distribution it should be 
possible to exert control over transverse and shear 
properties. | 

111) The wide variation in the literature of theoretical 
transverse Poisson's ratio certainly suggests an urgent 


need for accurate experimental results as well as a 


comprehensive survey of analytical techniques. 
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APPENDIX A: YIELD LOCUS CALCULATIONS 


The analytical method employed in this investigation 8 
requires that the macroscopic stress state of the composite 
be one of plane stress. It is also required that the rec- 
emaular form of the basic block, Figs. (2) and (3), be 
maintained Ns deformation. To ensure these results 
four distinct loading cases are superposed. These load 
cases and their boundary conditions are shown in Fig. (4). 
Load cases I and II are plane strain states with e 0. 


Load case III is a generalized plane strain case with 


E 


E 1.0. Load cases I, II, and III, are superposed to 


give Sy = 0. Load case IV is a state of longitudinal shear. 
All four cases have S xy = Syz = 0. Rectangular shapes of 
the basic block result from all load cases. 

Using finite element methods with linear strain triangles 
[2] the four loading cases are individually analyzed to 
obtaln the microstresses and their associated macrostresses. 
The microstresses are designated 914 with a subscript denoting 


the stress component and a superscript for the load case. 


For example, for load case I the microstresses obtained are 
I I I I 
σι, σαν, бу, ος (A-1) 
Similar results are obtained for problems II and III. The 


macrostresses for load case I are designated 


osa t ο 6 
шс Ен (А 2) 


and similarly for load cases II and III. 


S S 
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Load case IV yields microstress results for 


ту 
2 апа E 


and macrostresses 


IV 


IV 
эб апа Syz 


(A3) 


(n 


In the notation of Lin, Salinas and Ito [8] Ili dil roO 


stresses of a combined loading state may be related to the 


macrostresses as 


Kor = [A] {5} 
4х1 4х3 3х1 
апа 
ο) πι τ; 
2χ} 2x2 251 
Where 
T 
{o} = M ee 
T 
{5 } =< 545,5, > 
T 
{o} € Ox20y2? 
_ T ; 


M A352 A 


EU 33 


A41 A42 A44 
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(A-5) 


(A-6) 


(A-7) 


(A-8) 


(А-9) 


(А-10) 


`(А-11) 





[А] = (n) 


The matrix A may be related to the results of load 


cases I, II and III as 





ERU x A11 A12 413 
που тт тїт 
my oy су 22 23| |.I gl! git 
= OS x 
E TI ITLL 
пу ху “ху 329333 Е тет. 
S. S © А-13 
Y Y Y 
тт TILL 
A A 
ο О σ 
Ze 5 Ζ 42 43 gi 511 gill 
Zu 2 


For convenience this may be written 
[т] = [A] lv] (A-14) 
4x3 4х3 3x3 

The matrix A can then be evaluated as 


[a] = [т] +] (A-15) 


4x3 4x3 3x3 


Thus the combined microstresses, [c], due to any uniform 
transverse or longitudinal deformation may be evaluated by 


substituting the results of Eq. (A-15) into Eq. (A-5) as 


ο: | С) (A-16) 


4x1 4x3 3х3 3x1 


In a similar manner for load state IV and an analogous. 


load case V associated with 53, shear loading, the matrix 


A may be evaluated as 
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IV V 
O О 1 





ма х2 XZ IV 0 
[A] a Y y κα 
I L 
2x2 O yz Ove 0 SU (A-17) 


This may be written as 


К = πι BI" (A-18) 
2x2 2x2 2x2 ü 
From Eqs. (A-6) and (A-18) the microstresses, lo), may be 
related to any specified shear loading as 


{с} = [т] [+|-' {5} (A-19) 
2х1 2х2 .2x2 Eu 


The derivation of Eqs. (A-16) and (A-19) is independent 
of the load state. For the particular case of plane macro- 
scopic stress treated in this investigation the macroscopic 
stress vectors become 

T 


(S) = «S. O 


cre =" OS (A-20) 
Equations (A-16) and (A-19) may be evaluated for the point 
by point microstresses associated with a plane macroscopic 
load condition as specified by Eq. (A-20). 

Von Mises yield criterion is assumed to govern initial 


yield. This requires that 1 


Шо = О εσ ^-tGg. ^tf 


= = " A JE 2 m 
σχσν 9,0 7 ооё = Y (A-21) 
Y is the elastic limit of the material for the tensile 


stress state. 


sS 





Equations (A-16), (A-19), (A-20), and (A-21) may be 
combined to form a yield equation for plain macroscopic 
loading in terms of the macrostresses and the point depen- 


dent components of the A matrix. This equation is 


2 ΗΕ ο 2 2 _ _ _ 2 
С 2353508055299 02135237228^43 Άληλῃα) | 52 


an 
uc o3 521553 '^41^43? 


1 
gs (A4 1A23*A4 3331 35 A45*23h44 44/4 5*^43h))] Sz x 


4- 


2 l 2 2 zus = = 2 
hs σος is 221221 АА Ау 5х 


+ 


н 2 2 p ONE 
(A); А51) S ul Xo = 0 (A-22) 


For given values of 5 and $,, Eq. (A-22) may be solved 
for the value of S, necessary to cause yield. Each point 
within the fundamental block will give a different value of 
S, dependent on its microstress state. Each combination of 
amd S,, yields two values of S, at each point, a maximum 
and a minimum. The minimum of the maximum values and the 
maximum of the minimum values determine two yield points for 
Ener given values of S. and Ses Continuing in this manner 


the initial yield surface for a composite subjected to a 


macroscopic plane stress state is calculated. 
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APPENDIX B: CONSTITUTIVE EQUATIONS 
FOR A HOOKEAN FIBROUS COMPOSITE MATERIAL 


The well known equations of elasticity relating stress 
components and strain components in terms of two independent 
elastic constants represents a special case of generalized 
Hooke's Law for isotropic materials. Here we derive the 
Hookean stress-strain relations for a lamina oet anisotropic 
fiber reinforced composite material. 

The generalized Hookean stress-strain law relating stress 


and strain may be written in matrix form as 


Ex o Ju ls “ue li? 

Εν UE E975 8525 y 

Ez e a 25 

= (B-1) 

Exy 2119942 949 344 945. 9464. | Ps 

Ey z Em I 

€ a a a a a a S 

xz | 61 62 63 64 “65 66] | κα 


In a more compact notation this may be written 


le) = [a] {5 (B-2) 
6х1 6x6 6х1 


The Ej terms are engineering strain components and the Sij 
Stress components. 
Thus the generalized stress-strain relations contain 36 
unknown Hookean constants, 
ajj (у LAO 0) (B=3) 
Consideration of the strain-energy density function 27] 


establishes symmetry of the matrix of Hookean constants 
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and reduces the maximum number of independent elastic 
constants to 21. Further reductions in the number of 
elastic constants depends on the number of planes of 
symmetry present in the material. 

When the elastic constants at a point have the same 
values for every pair of co-ordinate systems which are 
Mirror images of each other in a certain plane that plane 
is called a plane of elastic symmetry. This means that if 
a certain plane, say the (x7 1X2) plane in co-ordinate 
system (X) X, X.) is a plane of elastic symmetry then the 


reflection of this system, (X) ,X5,X2), where 


2 
Um 
x) Xi 
кош =; 
ΙΗ (В-5) 
must have the same elastic constants. If you could actually 


"see" the elastic properties this means that the properties 
would appear identical at every point regardless of which 
side of the plane you were viewing. Alternatively Sines [28] 
suggests that if you could actually see these properties 
they would look the same even if the observer stood on his 
head, that is rotated 180° about his line of sight. This 
invariance is utilized to evaluate the elastic constants. 
Consider the reference system (X1 ,X5,X4) of Fig. (B-1). 
If the plane (1/2) is a plane of elastic symmetry a 
reflection through the х; axis may be carried out without 
effecting the elastic constants. This new reference system 


is designated by the primed axes of Fig. (B-1). 


J21 





Figure (B-1). 





x 


3 


XJ X1' 


Reflection of a Plane of Symmetry (X],X>). 


The general equations governing the transformation of 


Zx 


ΧΥ 


axes are 
a 2 2 2 
m = La u, e,+n] Eztk mje, M7 0] €y7 +0] %]€ 77% 
E. 2 2 2 
Ey = 23 *E¿Ho e, +n) Ez + L2mM)E tMmonge  „tna%ze 
E. 2 2 2 
g£" = ХЗ Ex +Mz ΕΥΠ Ez +L3M3€+M3N3Eyz+Hn343€ zx 
ο = Е € 4 TUNE 
Exy = 2°, 09°, +2mm2 "y t2njn5",t("4m4*m; 2) 
2l (m1n2?*njm5) €, ;* (nj * 5 * 105) €7. 
А 
Ενα = 22723€,+2m2mM3€,+2n203€7+ (22M3+M223) Ex y, 
+ (m5n3*n2m3) €, 7t (n5*4**5n3) €,. 
1 
Ezx = 2232] €x+2m3M7 €, 42030] €7+ (23M7+2]M3) Ex, 
+ (m3n,*n3mj) e, 7t (n321*£3n1) €zx 
for strain components and 
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(B-6) 





И — 2 2 2 
p - MSytmjSytnjS,t2*)miSy,*2mn1S, ,*2nj 5, 

B o2 2 2 
Ὃν = 225: +m28,+n757+247m257,+2m9098,, 74202425 zx 

Ш 2 2 2 
E = 235,+m35,/+n385 +2Ё#3ш35ху+2тзпз5,„+2п3Ё35уу 

T 

S (min5*njm2) $,,,* (n125*41n5) 87x (B-7) 
E = 22235,+m2m35,+m9n357+ (27m3+m243)S,,, 
+ (m3n3+n7m3)8,,,+(n223+%27n3)S,y% 

ENS = 


zx = 232157 +m3m1Sy+n3n757+(23m1,+m327]) Sy, 


+ (m3n,*n3mi)S,,*(n3g1t£3n1) $5. 
for stress components. The terms 
1:, mi, ni ИР) (B-8) 


are direction cosines defined by the following table: 





mmis indicates that 1, cos (x1,xj)etc. For the reflection 


of Fig. (B-1) the direction cosines are: 
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Using these direction cosines to evaluate the transformation 


relations Eqs. (B-6) and (B-7) yields 


— 1 ==. 
5х = 5х Ex 7 Ex 
= £.'-£ 
By Sy Y Y 
к 1 
52 = 5, Boc ασ 
1 1 
Sxy =Sxy “ху = в (B-9) 
n — = Ev) 
1 
Soy 7 zx бах = бах 


The stress-strain relation in the new co-ordinate system 


may be written 


T 
1 1 1 


T 
; ' 1 m | 
x 


The elastic constants are then evaluated by equating the 
transformed strains to the original strains in accordance 
with Eq. (B-9). Consider for example = πε In detail 


this requires that 

a118x+a125y+a3135772145xy"21558y2t9168%z = 
а115х+а125у+а135„+а145ху+а155у2+а1652х (B-11) 

The constants are invariant with transformation, requiring 


that 


-а1 4 Sxy = aja Sxy (B-12) 


For a non zero stress condition this can be true only if 
ais = 0 (B-13) 
Carrying this out for the other relations of Eq. (B-9) 


shows that 


а: -=а14=а24=а25=а. =аз5=а46=а-в=0 (B-14) 
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Thus a single plane of elastic symmetry reduces the 
number of independent elastic constants to 13. 

If a material has three orthogonal planes of symmetry 
a reflection of a second plane may be carried out without 
changing the Hookean coefficients. Reflecting the (X5,X3) 
plane through the Xi axis results in the transformed axes 


E Fig. (B-2). In this case the 





д 
NS 
En m 
хз ‚ХЗ t 
Figure (B-2). Reflection of a Plane of Symmetry (x2,%,) 


direction cosines are given by 





The transformation relations Eqs. (B-6) and (B-7) give 
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x x x x 
$ = с p c 
y Y Y Y 
' I 
= €E =€ 
S, S, И = 
Ü 1 
= = E E UE (B-15) 
ху ху ху ху 
S ' == E' = € 
y 2 yz yz yz 
©” = -S e" = HE 
ZX σος ZX 2% 


Using these relations it may be shown that 


ley tee 2745 7 0 ` (B-16) 


The matrix of elastic constants may now be written 


[А], = a33 0 0 0 (B-17) 
6x6 
244 0 0 
symmetric ass 0 
acc : 


where the subscript o indicates that this matrix is defined 
for an orthotropic material. There are now nine independent 
elastic constants. 

Further reductions in the number of independent elastic 
constants depend upon the number of planes of isotropy, if 
any, present in the material. A plane of isotropy is one 
in which the elastic constants are invariant with respect 
to direction within the plane. 

Assume that a plane of isotropy does exist, the (Χ1/Χ2) 


plane of Fig. (B-3). This indicates that the Hookean 
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coefficients are independent of the orientation of the 

(x, x.) axes within this plane. An infinite number of 
orientations could be chosen within this plane to assist in 
evaluating the remaining unknowns. A 45 degree rotation, 


shown in Fig. (B-3), was selected for ease of calculation. 





0 

ΧΩ ΧΑ 

Figure (B-3). 45 Degree Rotation of Plane (ху,х.) 
The direction cosines for such a case are: 


1 2 3 





Again using the transformation relations Eqs. (B-6) and 


(B-7) one obtains 


Ех Я Ч Ех + EV t 5 Ex, 
ES - Е = Εν - % Exy 
Es = Є, (B-18) 
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E | a 
Ey z == ү? Ey z a бух 
Е = Де E E 


S! = = З > sn 

Se = S, + š Suse 

D з, (B-18) 
Sxy = 5 Ux 7 А ὃν | 

Syz 5 L + خر‎ 5 2х 


' rl 1 
Sax ^ > Sy t Sox 


Using these relations it may be shown that 


a55 7 266 

A 222 

ауз = 453 (B-19) 
μα, ΠΗ, 


A material exhibiting one such plane of isotropy is 
called a transversely isotropic material. There are five 
remaining independent elastic constants and the matrix of 


elastic constants may be written 


ayj 312546 0 0 0 
E E о 0 
= a 0 0 -2 
fa " 33 0 (B-20) 
6x6 2(a,,7a3]2) 0 0 
symmetric ass 0 
a55 
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where the subscript TI signifies transversely isotropic. 
The stress-strain relations for this case may be written 


te} = [Apr] 15) 
6х1 6х6 6х1 


(B-21) 


One of the assumptions of this study is that a composite 
material may be considered to be homogeneous in a macro- 
scopic sense. This assumption allows the composite to be 
considered as a transversely isotropic material, and Egs. 
(B-21) are the constitutive equations for a composite 
material. Introducing the engineering constants these 


equations become 


ν 
Е = Vesey = Los 
x x T Ζ 
He Y Er 
1 Ур 
БЕ, S7) = 5 
y Em Y T x En Ζ 
Є = = VL Е А E S (B-22) 
2 E. X y E, Ζ 


uy cd 
_1 

yz Gr yz 
_1 
Sup 


LX 


En and Sy are the Young's modulus and shear modulus in 


the transverse direction, in the plane of isotropy. Er and 


G, are corresponding values in the longitudinal direction, 


perpendicular to the plane of isotropy. Poisson's ratio 
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Ут signifies the contraction in the plane of isotropy due 


to tension in the same plane. Poisson's ratio Vr, represents 
the contraction in the plane of isotropy due to tension in 


the longitudinal direction. The transverse shear modulus 


is a dependent quantity. 
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ας C: YIELD LOCI 


The figures included in this section represent the 
yield loci that were generated during the course of this 
investigation. Each plot is labelled with the constituent 
materials, filament volume fraction, Ver and EM and 5, 
scales. The outermost locus for all cases corresponds 
to Sxz = 0. The S,, increment between successive loci is 
also noted on each plot. 

Those plots that were used in the parametric analysis 
are labelled in terms of the parameter being investigated. 
All other properties are held constant. 

Note that in some cases the S, and S, scales are dif- 
ferent from one another. This occurs when the matrix 
materials are quite soft in comparison to the filament 
material. The transverse strength 1S so much less than 
the longitudinal strength that the surface plots nearly 
as a Straight line if the same scales are used for both 


axes. 
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Figure (C-25). 
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Figure (C-26). 
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Figure (C-27). 
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Figure (C-29). 





APPENDIX D: PROGRAM OVERLAY STRUCTURE 


The overlay root-segment structure is shown in Fig. 
(D-1) where each rectangle represents a subroutine. The 
subroutine functions are: 

MAIN -Remains in core and controls calling sequence 
during execution. 

BANDEC - Solves banded matrix problems. Used in ΜΑΙΝΙ 
and MAIN2. 

MULI -Multiplies a symmetric banded matrix with a 
vector. Used in MAINI1 and MAIN2. 

MAIN1 - Calculates microstresses and macrostresses for 
elastic loading of load cases I, II, and III. 

DAVES -Computes element stiffness and assembles system 
stiffness. 

CALTAU -Calculates element stresses and strains. 

REVGNK - Revises system stiffness to accommodate 
boundary conditions and applied loads. 

MAIN2 - Calculates microstresses and macrostresses for 
longitudinal shear elastic loading, load case IV. 

MAIN3 - Calculates Initial Yield Locus. 


MATINV - Matrix inversion with accompanying solution 


of linear equations. 
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